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SUMMARY 


Hydrodynamic measurements were made with a triaxial hot-wire in the 
full-coverage region and the recovery region following an array of injec- 
tion holes inclined downstream, at 30° to the surface. The data were 
taken under isothermal conditions at ambient temperature and pressure for 

two blowing ratios: M = 0.9 and M = 0.4. (M = p. U. /p U where U 

jet jet 00 00 

is the mean velocity and p is the density. Subscripts "jet” and M °°" 
stand for injectant and free stream, respectively.) Profiles of the three 
mean velocity components and the six Reynolds stresses were obtained at 
several spanwise positions at each of five locations down the test plate. 

High turbulent kinetic energy (TKE) levels were found for low blow- 
ing, with low TKE levels for high blowing. This observation is signifi- 
cant when coupled with the fact that the heat transfer coefficient is high 
for high blowing, and low for low blowing. This apparent paradox can be 
resolved by the hypothesis that entrainment of the main stream fluid in a 
3-D manner along the lanes between wall jets must be more important than 
the turbulent mixing in heat transfer for the high blowing ratios (close 
to unity) . 

The flow in the recovery region can be described approximately in 
terms of a two-layer model: an outer boundary layer and a two-dimensional 

(2-D) inner boundary layer. 

A one-equation model of turbulence (using TKE with an algebriac mix- 
ing length) was used with a 2-D computer program to predict the mean vel- 
ocity and TKE profiles in the recovery region. Mixing- length values 
calculated from the data were used in inputs to the program. The mixing- 
length distribution was assumed to be piecewise continuous, a heuristic 
fit to the data consistent with the concept of the two quasi-independent 
layers observed in the recovery region. This mixing length, along with 
a set of otherwise normal constants (for 2-D boundary layer predictions) , 
successfully predicted the spanwise-averaged features of the flow. 

A new hot-wire scheme was developed to make measurements in the three- 
dimensional (3-D) turbulent boundary layer over the full-coverage surface. 
The method uses a triaxial hot-wire and an analog device for real-time 


xi 



data reduction. It does not use the low-fluctuation assumption and can 
tolerate an unknown flow direction within limits (±30° for the mean 
velocity). The method is very fast and quite practical for taking large 
amounts of data in 3D, high-fluctuation turbulent flows. 
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Chapter 1 


INTRODUCTION 


1.1 Background 

Modern gas turbine technology requires high turbine inlet tempera- 
tures to increase thermodynamic efficiency, well above the usable temper- 
ature of the materials from which the turbine blades are manufactured. 

One way to protect the blades is by injecting a coolant through the blade 
surface, thus forming a cooling film between the hot combustion gases 
and the metal surface. There are several ways of injecting the coolant 
onto the blade surface. A comparison of different methods is given by 
Brunner (1969). Two of the more common are tra n spiration cooling through 
a porous plate surface and f ull-covei age film cooling through an array 
of closely spaced multiple rows of small holes, as discussed by Esgar 
(1971) . Transpiration cooling achieves lower heat transfer to the sur- 
face and has better aerodynamic performance than full-coverage film cool- 
ing but has problems with structural strength and is susceptible to 
plugging (Goldstein, 1971). Full-coverage film cooling seems more prac- 
tical, but not enough is known about the heat transfer and hydrodynamic 
mechanisms to allow accurate prediction. Such information is necessary 
before accurate designs can be developed. 

There have been many studies of heat transfer with film cooling. 

A general review is given by Goldstein (1971). A review of discrete-hole 
film cooling is given by Choe et al. (1975). Most of the early experi- 
mental research on full-coverage film cooling concentrated on measurement 
of film-cooling effectiveness with the heat transfer coefficients. The 
hydrodynamics was treated as secondary in importance. A review of the 
experimental work on this topic is also given by Crawford et al. (1976). 

Le Brocq , Launder & Priddin (1971) studied the effects on effective- 
ness, n , of hole-pattern arrangement, injection angle, ratio of densi- 
ties of the coolant and mainstream, and blowing ratio. Some profiles of 
mean velocity were taken near the holes, with a pitot probe. They found 
that the staggered hole pattern was more effective than the in-line 
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pattern. A critical blowing ratio was found to exist, which gave maximum 
effectiveness. Above this blowing ratio effectiveness decreased again. 
Slant-angle injection was found to be more effective than the normal in- 
j ection. 

Launder & York (1973) studied the effects of mainstream acceleration 
and turbulence level on film-cooling effectiveness, using a staggered 45° 
slant-hole test section. Again, only a small number of mean velocity 
profiles were taken. Acceleration seemed to increase effectiveness, but 
the free stream turbulence did not have a significant effect. 

Metzger, Takeuchi & Kuenstler (1973) studied both effectiveness and 
heat transfer on a full-coverage surface with normal holes spaced 4.8 
diameters apart, arranged in both in-line and staggered patterns. They 
did not make any measurement of hydrodynamics. Again, it was observed 
that a staggered pattern was more effective. 

Mayle & Camarata (1975) examined the effects of hole spacing and 
blowing ratio on heat transfer and film effectiveness for a staggered- 
hole array with compound-angle injection. It was found that for P/D = 

10 and 8, higher effectiveness was obtained than with P/D =14. No 
measurements were made of t lie hydrodynamics of the flow field. 

Choe et al. (1975) studied the effects on heat transfer of hole 
spacing, blowing ratio, mainstream velocity and conditions upstream of 
the discrete-hole array. Normal injection was used with a staggered 
array of holes spread 5 and 10 diameters apart. Mean velocity profiles 
were taken with pitot probe to obtain spanwise-averaged values. From 
the spanwise-averaged profile, mixing length was obtained and used in 
predictions of heat transfer data with a "zero-equation model" of turbu- 
lence. They observed that the P/D = 5 case with a comparable F was 
more effective than P/D = 10. Furthermore there was no pronounced 
effect in the initial blowing row. 

Crawford et al . (1976) also give a summary of analytical work in 

the field of full-coverage film cooling. Goldstein et al. (1969) and 
Ericksen, Eckert & Goldstein (1971) used superposition of film-effective- 
ness data for individual jets to predict n by modeling injection as a 
point heat source. Mayle & Camarata (1975) developed an improved super- 
position method to predict their full-coverage data. Pai & Whitelaw 
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(1971) and Patankar, Rastogi & Whitelaw (1973) investigated the predic- 
tion of wall temperature and effectiveness downstream of 2- and 3-D film- 
cooling slots. For the 2-D slot injection, the boundary layer equations 
were used along with an augmented mixing length model to represent the 
effect of injection. For the 3-D injection, 3-D "Navier-Stokes equations 
were solved numerically by reducing them to elliptic form in the lateral 
plane and to parabolic form in the streamwise direction. 

Herring (1975) used a finite-difference method for predicting the 
flow over a full-coverage film-cooled surface. Lateral averaging in the 
full-coverage region was invoked to justify using 2-D boundary layer 
equations. Terms arising from the spanwise averaging were obtained from 
consideration of the jet-boundary layer interaction. Predicted velocity 
profiles were reported but showed some problems in the initial regions 
of injection near the wal3. . 

Choe et al. (1975) developed a finite-difference method for predict- 
ing heat transfer with full-coverage film cooling, solving the two- 
dimensional boundary layer equations. (These equations have a form simi- 
lar to those given by Herring, 1975.) Choe et al. (1975) used the concep 
of local averaging, with a different model for the injection process, 
the nonlinear terms, and augmented turbulent mixing. By using an augmen- 
ted mixing length in a zero-equation model for turbulence, he success- 
fully predicted most of their Stanton number data for low and moderate 
blowing in the full-coverage region. The predictions in the recovery 
region and for high blowing were less accurate. 

Crawford et al. (1976) repeated the conditions covered by Choe et 
al. (1975); however, they used a different full-coverage surface: 30° 
slant-hole injection with a staggered array of holes, with two different 
hole spacings (5 and 10 diameters). For each injection, Stanton number 
data were taken for two values of injectant temperature, corresponding 
to 6 equal to zero and one. Linear superposition was then used to ob- 
tain Stanton numbers as a continuous function of the injectant tempera- 
ture. Stanton numbers were measured for a range of injection mass flux 

5 

ratios (M from 0 to 1.3) and Reynolds numbers (Re^ from 1.5 x 10 
to 5 x 10 6 ). 
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With the injection temperature equal to the wall temperature, the 
Stanton number decreased below the normal flat-plate value and reached a 
minimum at M of 0.4. A higher M caused an increase in the Stanton 
number. For the recovery region, downstream of the five-diameter hole 
array, two distinct data trends were observed. For low M the Stanton 
number began to recover immediately from the effects of blowing, while 
for high M the Stanton number either remained constant or dropped 
throughout the recovery region. This latter behavior suggested that 
interrupted hole patterns may be an effective cooling scheme. 

In Crawford’s work (Crawford et al., 1976) both integral and dif- 
ferential analyses were carried out on the data sets. A differential 
prediction method was developed which used a two-dimensional boundary 
layer program. The injection process was modeled as slot injection, and 
the turbulence was modeled with a zero-equation model, using an augmented 
mixing length obtained from the spanwise-averaged mean velocity profiles. 
Mean velocity profiles were obtained with a pitot probe. Downstream of 
the jets the profiles may be in error because of the static-pressure 
islands created around the jets. In Crawford’s data, most of the five- 
hole diameter data were successfully predicted in the full-coverage re- 
gion. The predictions in the recovery region and for high blowing ratios 
were not very good, as was also the case with Choe et al. (1975) . 

Very little of the literature mentioned above gives any information 
about flow field measurements and, in particular, measurements of turbu- 
lence quantities. In the following papers, however, some attempts were 
made to make measurements of hydrodynamics. Ramsey & Goldstein (1971) 
measured the mean velocity and turbulence intensity profiles ( no t the 
turbulent kinetic energy) for normal injection downstream of a heated jet 
at blowing ratios of 1.0 and 2.0. Metzger. Carper & Warren (1972) 
measured the mean velocity profiles both upstream and downstream of a 
two-dimensional injection slot. Hartnett, Birkebak & Eckert (1961) mea- 
sured boundary-layer velocity profiles for a number of positions down- 
stream of injection through a single slot. Foster & Haji-Sheikh (1974) 
measured mean velocities downstream of flush, normal-injection slots. 
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I 

In these references it was frequently observed that the literature 
is quite weak from the point of view of hydrodynamic data, especially 
for full-coverage surfaces. 

1.2 Objectives 

This research represents an extension of the work by Crawford et al. 
(1976), whose experiments consisted mainly of heat-transfer coefficient 
measurements. As explained, the predictions of Crawford were good for 
low and moderate blowing ratios in the full-coverage region, but problem- 
atic in the recovery region and for high blowing ratios . It was felt 
that the key to understanding these result lay in the structure of the 
turbulence. 

A decision was made to adapt a one-equation model of turbulence for 
prediction of the hydrodynamics within the recovery region, since the 
zero-equation model did not give satisfactory results. It was thought 
that a higher-order model might better represent the real process. 

The objectives of this research were twofold: 

• to make 3-D measurements of the turbulence structure and mean 
velocity field in the 3— D turbulent boundary layer over the full- 
coverage film-cooled surface and in its recovery region, 

• to apply a one-equation model of turbulence for predictions of the 
recovery region hydrodynamics with the empirical input supplied by 
the experiments. 

1. 3 Outline of the Work 

A new hot-wire anemometer method of measurement was devised to allow 
measurements of the turbulence structure in the 3-D turbulent boundary 
layer over the full-coverage surface. Due to the interaction between the 
mainstream and the jets, a complex flow field is created in the full- 
coverage region where the local flow is 3-D in the vicinity of the jets, 
the local flow direction is unknown, and the fluctuations are high. 

The new hot-wire method is believed to be an important contribution. 
The present research could hardly have been done without it. The method 
is described in Chapter 3 of this thesis. 
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The physics of the recovery region flow are described, in Chapter 4, 
as well as the predictions of the recovery region hydrodynamics using 
the one-equation model of turbulence. The physical modeling is accom- 
plished by modeling the length scale of turbulence. The length-scale 
model introduced in this chapter was successful in predicting the flow 
with the one-equation model. 

Chapter 5 presents the data taken on the full-coverage surface and 
discussions of the data. The volume of the data taken was so large that 
only typical results are discussed. As far as the author knows, these 
data comprise the first detailed study of this process. There will be 
much to explore in the data beyond the scope of the present report. 

Chapter 6 contains the summary and recommendations about the com- 
plete research. 

The building of the 3-D Turbulent Flow Analyzer was funded by the 
Thermosciences Affiliates. The Scientific and Technical Research Coun- 
cil of Turkey (Turkiye Bilimsel ve Taknik Arastirma Kurumu) provided all 
the first author’s living expenses and tuition during his four-year stay 
at Stanford University. 
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Chapter 2 


APPARATUS AND INSTRUMENTATION 

The Discrete Hole Rig is the hear transfer apparatus which carries 
the full- cover age film-cooling test section. It is a closed loop wind 
tunnel which delivers air at ambient pressure and constant temperature. 
This system was previously used by Choe et al. (1975) and Crawford et 
al. (1976) . It is described in detail in their work. Only those aspects 
relevant to the present research will be given here. A flow diagram of 
the systems which make up the rig is shown in Fig. 2.1; its photograph 
is given in Fig. 2.2. 

2 . 1 The Test Surface 

The test surface of the Discrete Hole Rig is made from copper plates. 
It has three parts: a preplate, a full-coverage region, and an after- 
plate. 

The preplate is made of 24 individual copper plates, each about 
2.6 cm long in the flow direction and 46 cm wide. 

The full-coverage section is composed of 12 copper plates, each 0.6 
cm thick, 46 cm wide, and 5 cm long in the flow direction. The first 
plate does not have any holes, but the remaining 11 downstream plates 
plates have alternately nine and eight holes each. The holes are each 
1.03 cm in diameter and are spaced on five-diameter centers to form a 
staggered hole array. A photograph of the full-coverage section can be 
seen in Fig. 2.3. 

Delivery tubes for secondary air extend beneath the surface at a 30° 
angle to the plate surface, giving a 30° slant -angle injection in the 
direction of the main flow. A schematic of the full— coverage surface can 
be seen in Fig. 2.4. 

The afterplate is identical to the preplate. This is the recovery 
region for the flow after the blown part of the test section. See Craw- 
ford et al. (1976) for more detailed information about the test surface. 
Twenty-four instrumented strips in the afterplate permit measurement of 
the Stanton number distribution in the recovery region. 
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2.2 Instrumentation 


This section will give information on the instruments used for mea- 
surements in the turbulent flow field. The amount of detail given about 
a certain instrument depends on the importance of its use in the measure- 
ments. Detailed information is given about the calibration of the tri- 
axial hot-wire probe. 

2.2.1 Temperature 

Since the experiments were performed under isothermal conditions and 
the purpose was to make measurements of hydrodynamics, there was no need 
for extensive temperature measurements. Only the free-stream temperature 
was measured. This measurement was made with an iron-cons tantan thermo- 
couple whose junction was parallel to the flow. The thermocouple probe 
did not require any special shielding or precautions because the operat- 
ing temperatures were about ambient and the velocities were low. An in- 
tegrating digital voltmeter (HP 2401C) was used to read out the thermo- 
couple output. 

2.2.2 Pressure 

Tunnel static pressure was measured with an inclined Meriam manom- 
eter using a 0.824 specific gravity fluid. The manometer had a range of 
0.5 in with a smallest division of 0.005 inches of water. Static pres- 
sure was measured from the taps located in one of the tunnel sidewalls. 
Their locations and properties are given by Choe et al. (1976). The main- 
stream dynamic pressure was measured with a kiel probe corrected to a 
Combist micromanometer whose smallest division was 0.0005 inches of manom- 
eter fluid of specific gravity 0.82. The tunnel geometry was adjusted to 
achieve a uniform static pressure throughout the test section. 

2.2.3 Secondary Air Flow Rate 

Flow rates through the secondary air delivery tubes were measured 
with hot-wire flow meters permanently built into the rig. More informa- 
tion about the flow meters is given by Choe et al. (1975). The digital 
voltmeter mentioned in section 2.2.1 was used to read out the flow meter 
outputs. Calibration of the meters w 7 as done individually, using reference 
flow meters. 
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2.2.4 Turbulent Flow Field 


a. Instrumentation 

The main instrumentation for the experiments was the hot-wire system 
used in the measurements of the turbulent flow field. It consists of a 
DISA claw- type triaxial hot-wire probe, three constant temperature ane- 
mometers (TSI 1050), three polynomial linearizers (TSI 1052), and the 3-D 
Turbulent Flow Analyzer. 

An HP-2100 digital computer was used for data acquisition and reduc- 
tion. 

An integrating digital voltmeter (IDVM) (HP-2401C) and an RMS meter 
(TSI 1076) were used occasionall}' to reduce and read the outputs of the 
3-D Turbulent Flow Analyzer instead of the computer. 

More detailed information about the hot-wire instrumentation is 
given in Chapter 3. 

b . Calibration of the triaxial hot-wire probe 

The triaxial hot-wire probe was calibrated on the calibrator de- 
scribed by Pimenta et al. (1975). This calibrator produced a uniform flow 
of air at constant but adjustable temperature, with moderate- to-low turbu- 
lence levels (about 0.8%). 

The probe was calibrated in the free jet at the exit of the calibra- 
tor. After the flow in the calibrator was adjusted to the desired tem- 
perature and the anemometers were set, the following procedure was applied 
for the process of calibration. 

The axis of the probe was aligned with the flow direction by trial 
and error. To check the alignment, the output of each anemometer was re- 
corded at a certain orientation of the probe. Then the probe tip was 
rotated 180° (the probe holder was designed so that the probe tip could 
be rotated to 12 positions, each 30° apart) and the anemometer outputs 
were recorded again. If the probe was aligned properly, the outputs of 
the same anemometer at 0° and 180° positions of the probe tip should have 
read the same (for 10 seconds integration of the output signal, the ac- 
ceptable value is x.xxx ± 0.002 volts for a velocity about 10 m/sec). 

If the numbers were not the same, the probe position was readjusted until 
the alignment criterion was satisfied. It usually took about five trials 
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to obtain the correct alignment. The final alignment position was within 
1-2° of the first visual alignment. During this procedure it was very 
important to make sure that any change in readings was not caused by a 
change in either the temperature or the flow rate. 

After the probe was aligned with the flow direction the calibration 
started. First the flow speed was adjusted to the desired value by ad- 
justing the valve of the calibrator. The flow speed was measured 
through a static pressure tap located on the plenum chamber where the 
pressure is equal to the total pressure of the jet at the nozzle exit. 

The pressure was read with the Combist micromanometer mentioned in sec- 
tion 2.2.2. Afterwards, the output of the three anemometers was read 
with the IDVM, described before. Ten seconds of integration of the out- 
put was enough for the calibration. 

It was necessary to check the temperature during calibration. This 
was done by checking the cold resistance of one of the wires from time to 
time. 

Finally, the obtained values of static pressure and anemometer out- 
puts along with the temperature, ambient pressure, and humidity were fed 
to a computer program prepared by Coleman (Coleman et al. , 1975) to ob- 
tain the linearized calibration equation. 

The typical values used for some of the variables during the cali- 
bration were the following: 

Cold wire resistances: - 9.5 ft 

Overheat ratio: - 1.5 

Maximum change in the cold resistance 

during one calibration: - ± .02 ft 

Temperature range (for different cali- 
brations): 18°C-26°C 

Velocity range during one calibration: 4-22 m/sec. 
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Figure 2.1 Flow schematic of wind tunnel facility, the Discrete Hole 
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Cross-sectional drawing of the full-coverage section. 



Chapter 3 


REAL-TIME HOT-WIRE MEASUREMENTS IN THREE-DIMENSIONAL FLOWS 

Hot-wire anemometers are the most widely used devices for making 
measurements in turbulent flow fields. Satisfactory results are obtain- 
able in one- and two-dimensional flows whose directions are known and 
whose turbulence structures are stationary and of low amplitude (typi- 
cally on the order of 10% of the mean velocity). Three-dimensional flows 
of unknown flow direction, with high fluctuations or large intermittent 
structures, cannot be measured accurately with classical techniques, and 
even the stationary, low-fluctuation cases are very tedious. 

The first step in most classical methods involves finding the time- 
averaged value of the fluctuating component of the anemometer output. 
Details of the turbulence structure are deduced by solving sets of simul- 
taneous equations using measurements at different angular positions. 

These equations are formed by expanding the expression for effective 
velocity in a Taylor's series around the mean and truncating at second 
order (usually) . When the equations are weakly convergent (as when many 
angular positions are used), the uncertainty in the calculated values of 
the turbulence quantities is high, especially for the higher-ordered 
terms . 

The object of this work was to develop a real-time method for hot- 
wire anemometry which would be applicable to three-dimensional flows of 
unknown flow direction (within limits) and with high fluctuation levels. 

3 . 1 B ackgro un d in Hot-Wire Measurements 

3.1.1 Cla s sical ti me-ave raged reduction s 

Most of the existing hot-wire methods are good in 2-D flows with low 
fluctuations and known flow direction. Almost all use either single 
horizontal wires, single rotatable slant wires or cross wires. The usual 
practice in data-taking is to align the probe with the known flow direc- 
tion, to read the mean output of the anemometer or linearizer with a 
digital voltmeter, and the average fluctuation value with an RMS meter. 
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The desired quantities, mean velocities and Reynolds stresses are deduced 
from the truncated Taylor series expansions of these effective quantities. 
The truncation and solution process requires some simplifying assumptions, 
in particular, that the fluctuations are small compared to the mean values. 

When the flow direction is unknown, the data-reduction process is 
quite complicated. To explain this point better, the relation between 
the effective velocity and the velocity components is given in Eqn. (3.1) 
taken f rom Jorgensen (1971) , with the nomenclature shown in the following 
f igure . 



2 2 2 2 2 2 
U., = X + k n Y + k Z (3.1) 

err i z 

Here, X, Y, Z are the velocity components in the directions of the wire 
coordinates x', y T , and z', respectively, and k^ and k 2 are pitch 
and yaw sensitivities of the wire. The following figure shows a single, 
rotatable slant wire and the wire orientation with respect to the labora- 
tory coordinates x, y, and z. 
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In this figure <f> is the slant angle, in the plane of prongs, between 
the wire and the y' direction. The angle <j> is zero when the wire is 
perpendicular to the axis of the prongs. The angle 0 is the roll angle 
in the (yz) plane, and is zero when the prongs are in the (xy) plane. 
Velocity components U, V, and W are in the laboratory coordinates, in 
the x, y, and z directions, respectively. With this nomenclature, 

Eqn. (3.1) can be written in terms of velocities in laboratory coordinates, 
as follows: 


U 2 = A(U 2 ) + B(v 2 ) + C(W 2 ) + D (UV) + E(VW) + F(UW) (3.2) 

err 


where 


A = cos 


2 2 
+ k^ sin 


2 2 2 2 2 2 

B = (sin 4> + k 1 cos <f>)cos 0 + sin 0, 

2 2 2 2 2 2 

C = (sin <f> + k^ cos (f>)sin 0 + k^ cos 0, 


(3.3) 


D = (l-k^)sin 2<f> cos 0, 

2 2 2 2 

E = (sin <J) + k^ cos <j> - k^)sin 20, 

2 

F = (l-k^)sin 2<f) sin 0. 

In Eqn. (3.2), each velocity can be written as the sum of the mean and 
fluctuating parts: 

U = U + u* 

V = V + v’ (3.4) 

W = W + w’ 

and 

U „ = U __ + u' (3.3) 

eff eff eff 


In the classical time-averaged method, the quantities recorded are 

— Y~ 

U and u’ _ for each of several angular positions. If the first- 

eff eff 

order fluctuation descriptions are substituted into Eqn. (3.2), a very 
complicated equations results which can be brought to a theoretically 
solvable form by a Taylor series expansion of both sides. This expansion 


17 



is valid only when the flow has a strongly preferred direction and low 
fluctuations (i.e., V, W, and u’, v', and w T are each at least one 
order of magnitude smaller than U) . The resulting equations are given 
below: 


U = — - — f(2AU 2 + BV 2 + CW 2 ) + (DUV + EVW + FUW) 

ef f 2/K U L 


(3-6) 


,2 


+ (Au' + By’ + Cw’“) + (DuV + Ev T w ' + Fu'w T ) 


+ 0(3) 


& 


4A 


4A 2 + t ^ +f 2^ + 44 B V +4jF ^ + 2DF ^u ' 2 

u 2 u 2 u u u 2 

—2 
2 V 


+ {D 2 + 4B 2 — - + E 2 "■ + 4BD z + 2DE ^ + 4BE ^ ) v T 2 
U U U U U ' 

+ IF 2 + E 2 ~ + 4C 2 ~ + 2EF Z + 4CF Z + 4CE zz)w ’ 2 

u u u u u 

— 2—2 — — 

+ ( 4 AD + 4BD ~ + 2EF ~ + ( 8 AB + 2D 2 ) ~ + (4AE + 2DF) ~ 

u 2 U U U 


(3.7) 


^4 A 


V W 


_9 

V* 


w 


2 


V 


+ (2DE + 4BF) Z3-) u’v’ + 4AF + 2DE + 4CF + ( 4AE + 2DF) £ 

U ' \ u 2 u 2 U 


+ (SAC + 2F 2 ) ~ + (2EF + 4 CD) u’w’ + ^2DF + 4BE ~ + -CE — 
U U ’ ' U u 


Z 2 

I 


+ (2DE 4- 4BF) ~ + (4 CD + 2EF) ^ + ( 8 BC + 2E 2 ) 


2/ 


v w 


+ 0(3) 


Even with these assumptions the resulting equations are quite com- 
plicated. In two-dimensional flows with low fluctuations, where the flow 
direction is known, aligning the probe with the flow direction sets V 
and W to zero, allowing the simplified equation to be used. In cases 
where there is symmetry (u’w’ = 0 and v’w T = 0 ), the equations reduce 
to the following form, which is generally used for rotatable slant wires: 


U ff = SR U + 0(2) 


(3.8) 


,2 


,2 D“ ,2 F f 

1 4- \T 1 -i_ T.7 ’ 


u‘ iff = A u + 4 X v + 4 A w? + D uTv ’ + °( 3) (3.9) 
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One measures and u ^ff i- n eac h of four different positions, thus 

obtaining enough equations to solve for these unknowns. 

There are several ways to make measurements in 3-D turbulent flows 
using the classical time-averaged equations. The most important problem 
in 3-D flows arises when the mean flow direction is unknown. This direc- 
tion can be found with one of the following ways. Johnston (1970) mea- 
sured the local pitch and yaw angles for the mean velocity with a Conrad 
probe and then measured the Reynolds stresses with a horizontal wire and 
a rotatable cross wire aligned with the flow. The low-fluctuation assump- 
tion was invoked. Moussa & Eskinazi (1975) tried to measure the mean 
flow direction by using a rotatable slant wire. Making use of the direc- 
tional properties of hot wires, they calibrated the probe for all possible 
angles and prepared detailed charts which included the flow angles as 
functions of four mean voltages obtained at different rotations of the 
probe. Delleur (1966) used a cross wire to measure the flow direction, 
arguing that the cross-wire technique was twice as accurate as the single- 
wire technique. Both methods require the use of a hot-wire calibration 
curve for flow-direction measurement and consequently require frequent 
calibration of the hot wire to renew the calibration charts. 

Some other 3-D hot-wire methods were developed which do not require 
the flow direction to be known. Mojola (1974) gives a hot-wire method 
(rotatable slant wire or cross wire) for measuring the three mean compo- 
nents of the velocity and six Reynolds stresses without knowing the flow 
direction. His equations are valid when the probe is approximately 
aligned with the flow (i.e., a strongly preferred mean flow direction) 
and has only low fluctuations. Hoffmeister (1972) describes a scheme 
which employs a single rotatable slant wire to obtain three mean veloci- 
ties and six Reynolds stresses. In this scheme the interpretation of 
the anemometer voltages is based on calibrations of the probe over the 
entire range of angles between the wire and the flow which may exist dur- 
ing measurements. 

None of the preceding methods is practical for taking large amounts 
of data. Further, the accuracy with which the higher-ordered terms can 
be measured is seriously limited. As the number of terms retained in 
Eqn. (3.7) increases, the number of independent realizations required 


IS 



increases; thus the number of probe rotational positions increases and 
the strength with which the equations converge diminishes. 

The problem lies, basically, in the time-averaged approach to tur- 
bulence measurement. With one or two wires, one does not have enough 
information to solve for the instantaneous velocities, and, hence, time- 
averaging is required. With three wires, however, the whole problem is 
resolved and one can deal directly with the instantaneous velocities. 

3.1.2 Real-time reduction 

A real-time data-reduction scheme was used by Zimmerman & Abbott 
(1975). Using a triaxial hot-wire probe and an analog device, they 
solved for the three unknowns (three instantaneous velocities) in the 
effective velocity expression, in real time. The outputs of their analog 
device were the three mean velocities, the three instantaneous fluctua- 
tions, and the means of three correlations (off-diagonal components of 
Reynolds stress tensor) . This scheme was used to measure the Reynolds 
stresses in a 3-D turbulent boundary layer created by a skewed-leading- 
edged flat plate in zero pressure gradient air flow, which yields high 
cross-stream gradients of the turbulent quantities. 

3 . 2 The Present Approach 

The present approach also uses the real-time data-reduction idea. 

In fact, this work was finished at the same time as Abbott & Zimmerman's 
(1975) work, but the publishing of results has been delayed due to un- 
foreseen circumstances. 

Equation (3.1) introduced the effective velocity expression. The 
three unknowns are the three instantaneous velocity components in wire 
coordinates. To solve for these velocities in real time, one needs three 
simultaneous and independent values of the U ^ at one P°i nt in the 
flow field. Probes are available which have an orthogonal array of wires, 
such that the x* axis for one is the y* and z 1 axis for the others. 
The probe used in the present studies was a DISA triaxial-wire probe. 

The DISA triaxial wire probe has three wires each of which, has two 
leads and each of which is driven by a separate anemometer. The wires 
are mutually orthogonal, forming a right-angled coordinate system. The 
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sensors provided by DISA are 3.2 mm long, with 5 micron Pt-plated tungsten 
wires having copper- and gold-plated ends leaving a sensitive length of 
1.25 mm. In the present work, bare Pt-plated tungsten wires were used 
without gold plating, giving 3.2 mm active length. The three wires form 
a cone of apex angle 70.6° around the axis of the probe stem. The wire 
coordinate system can be seen in Fig. 3.1. With the special prong struc- 
ture of the probe, the effective velocity indicated by each wire is rela- 
ted to the velocity components in the wire coordinates in the following 
manner. 


U 

U 

U 


2 

ef f^ 

2 2 2 2 
k ll X + Y + k 2 l Z 

2 

ef f 2 

2 2 2 2 
k 22 X + k 12 Y .+ Z 

2 

ef f 3 

X 2 + k 23 Y 2 + k 2 3 Z 


2 

2 

2 


(3.10) 


and the linearized effective voltages (linear izer outputs) are related to 
the effective velocities as follows: 
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eff l = 

A 1 
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B l E eff 


A n 

+ 

B„E 

ef ^2 
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2 eff 


A„ 

+ 

B E _ _ 

ef f 3 

3 


3 eff 


1 

2 

3 


(3.11) 


In Eqns. (3.11), the A’s and B’s are constants obtained from calibra- 
tions of the wires. 

Equations (3.10) have three unknowns — the instantaneous velocities 
in the wire coordinates — which can be obtained from the equations shown 
below: 
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(3.12) 
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(3.13) 


M 


Once the instantaneous velocities in the wire coordinates are obtained, 
the instantaneous velocities in the laboratory coordinates can be ob- 
tained easily with a transformation of coordinates. 


"u“ 


"x" 

V 

= N 

Y 

_w_ 


_Z_ 


(3.14) 


where N is the coordinate transformation matrix from wire coordinates 
to laboratory coordinates. With the three instantaneous velocity compo- 
nents available in laboratory coordinates, the mean values and the Rey- 
nolds stresses can be obtained by using digital voltmeters, RMS meters 
and turbulence correlators, or with a digital computer having a fast 
enough A/D (analog-to-digital) converter. 

3 . 3 The Three-Dimensional Turbulent Flow Analyzer 

An analog device was built to solve the equations (3.13) and (3.14) 

using high-speed analog components. All quantities were magnitude-scaled 

using (U , the value which resulted in the maximum output of the 

° eff max 

linearizer. (For the present case (U ..,-) can be obtained by substi- 

r eff max 

tuting 10 volts for E i.e., (U ff ) - A + B x 10.) Also, in 

6 eff eff max 

calibration the problem is scaled such that the output in the real experi- 
ment will be between 4 and 8 volts, never reaching the limits of 0 volts 
and 10 volts. Under these normalization conditions it can be shown that 
none of the quantities in the equations to be solved goes over 10 volts 
during an experiment. After normalization, the equations to be solved 
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by the analog system are (3.13) and (3.14), with all the quantities nor- 
malized on (U j-j.) 

eff max 

A and B values can be set to be the same for all wires during 
calibration, since they are functions only of the minimum and maximum 
velocities of the calibration range and the voltages desired. If the 
stem of the three-wire probe is aligned with the calibration flow, all 
the wires are exposed to the same minimum and maximum velocities. 

The coefficient matrix M can easily be obtained using the known 
values of and k^ for each wire. For example, for k^ = 0.15 and 

1^2 = 1.02, the matrix is: 


M = 


~-0.5058 

0.5155 

0.4750 


0.4750 

-0.5058 

0.5155 

(3.15) 

_ 0.5155 

0.4750 

-0 . 5058_ 


used for 

k^ and 

k 2 are the same for all 


three wires. Uniform values of k^ and k£ were used in the present 


study, but the analog system is designed such that it can easily handle 
different values for k^ and k£ for each wire. An uncertainty analy- 
sis showed that a 1% change in k^ causes a change of only 0.02% in the 
indicated velocity, while a 1% change in k^ yields a change of 1.04%, 
regardless of the flow angle. Thus the final results depend only slightly 
on k^ but nearly linearly on the value of k^- As was shown by Jorgen- 
sen (1971), the change in k£ with flow angle is very small. 


N = 


cos GO 
0 

sin go 


-sxn a) 
0 

COS GO 


0 


cos a 


-sxn a 


0 

sin a 
cos a 


“ /3/3 /3/3 /3/3~ | 

-/6/6 -/6/6 /6/3 

-/2/2 /2/2 0 


(3.16) 


where the angles a and go are shown in the following figure: a is 
the angle of rotation around the probe axis and go is the angle between 
the probe axis and the x axis when the probe body is in the x-y plane 
The angles are measured in the direction indicated by the arrows. Angle 
a is taken to be zero when the third wire is in the x-y plane and the 
longer prong is below the shorter one. 
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The 3-D Turbulent Flow Analyzer has been designed and built to take 
the linearizer outputs E e ^ , E e ff^, and E e ff^ as inputs and give 
instantaneous (J, V, and W (the normalized velocities in laboratory 
coordinates) as outputs. Fig. 3.2 shows a photograph of the two control 
panels of the 3-D Turbulent Flow Analyzer. On the lower panel, there is 
one potentiometer for setting the A coefficient, and three potenti- 
ometers for setting the B coefficients. Three potentiometers are needed 
for B values because of circuit requirements. The nine coefficient po- 
tentiometers for the inverse Jorgensen Matrix (M) are shown, each with 
its own pair of test points for checking values. The meter provides a 
continuous display of R/C, the time-averaged value of the magnitude of 
the normalized velocity vector. The outputs are the normalized velocity 
components in wire coordinates X, Y, Z and a value of R/C, the root 
sum square of the components, divided by three. The upper panel accepts 
the X, Y, Z as inputs and contains nine potentiometers whose values can 
be set according to the wire position to calculate the values of U, V, 
and W, the normalized velocity components in laboratory coordinates. 

A flow diagram of the complete system is shown in Fig. 3.3. 

The electrical performance of the 3-D Turbulent Flow Analyzer was 
checked up to 20 kHz for phase shift and attenuation of magnitude. 
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Maximum phase shift was between 1 and 2°, and maximum attenuation of 
signal magnitude was approximately 0.1%. 

3 . 4 Processing of the Data 

There are several options which may be exercised once the instanta- 
neous values of the three components of velocity are available as analog 
electrical signals. Analog multipliers could be used to form instanta- 
neous products, with integrating voltmeters for time averages and RMS 
meters for fluctuation measurements. In the work reported here, the three 
signals were digitized and processed through an HP-2100 system. The A/D 
converter used had a capability of 45,000 points per second with which to 
process all three channels. This would have limited the acceptance to 
about 7500 Hz, using two points per cycle. Peripheral limitations on the 
current system in fact limited the response to 5882 Hz. Thus the high- 
frequency events cannot be seen through the present system, but the re- 

/ 72 

striction is in the A/D system, not the analog system. Values of u 
measured by analog and digital methods agreed to within 2-3%; thus the 
frequency limitation on the digital system seems not to have been impor- 
tant. Data were recorded onto digital tape and processes in computer- 
idle time, overnight. Twenty-two seconds of data-taking per point re- 
quired 11 minutes of post-acquisition processing on the HP-2100 to produce 
the values of the three mean components and the six Reynolds stresses. 

With all-analog processing, the nine outputs could be generated in what- 
ever time was felt necessary to achieve stationary values, based on the 
flow. 

The functions of the Flow Analyzer could all be performed digitally, 
by passing the three linearizer outputs through A/D converters by exe- 
cuting the matrix-inversion processes digitally. Using the same HP-2100 
system, with the same A/D conversion unit, the total time required for 
all-digital processing was 75 minutes to process 22.6 seconds of data- 
taking. The use of analog operations for the two matrix inversions re- 
duced the computing time by a factor of 8.5. Different computer systems, 
with different A/D conversion and higher speed would surely make the di- 
rect digital system look more attractive than these results. Whatever 
digital system is available, however, the two analog stages presently 
used will continue to be attractive from a cost-of-running standpoint. 
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3.5 Problems Created by the Size of the Triaxial Wire Probe 

The triaxial wire probe creates some problems in the measurement 
because of its relatively large size compared to the single- and cross- 
wire probes. The probe has three wires, with their centers located on, a 
sphere of diameter 3 mm to prevent interference due to the wakes of wires. 
Measurements made on this sphere must be attributed to the center of the 
sphere. When making measurements in high shear flows, there may be some 
uncertainty about where the effective measurement center is, since it 
will change depending on the orientation of the wires with respect- to the 
gradient. Another problem is the overall size of the probe, about 0.8 cm 
in diameter. It is not possible to obtain data closer than 0.4 cm to the 
wall. 

3 . 6 Qualification Tests 

The Thermosciences Laboratory in the Mechanical Engineering Depart- 
ment of Stanford University has a two-dimensional channel which gives 
fully developed mean-velocity and turbulence profiles at least to the 
second-order turbulence quantities. The performance characteristics of 
this channel have been fully explored; it has been used by several re- 
cent experimenters to calibrate their hot-wire technique, their probes, 
and their systems. The Flow Analyzer was tested in this channel to qual- 
ify its performance and to explore its limitations for turbulence measure- 
ments . 

For the qualification tests, the probe was mounted in a two-axis 
probe holder so that it could be rotated around the streamwise axis (roll 
angle a) and also tilted (pitch angle oj) against the approaching flow, 
as well as traversed to several different distances from the wall. By 
measuring at several distances from the wall, the system performance was 
recorded both for high shear regions (near the wall) and zero shear re- 
gions (at the centerline of the channel). The outputs were compared with 
the outputs of the other acceptable methods of measurement in the channel 
(single horizontal wire, pitot tube measurements, and linear shear stress 
distribution calculated from the pressure gradient along the channel). 
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The two-dimensional channel is 6.35 cm wide and 117 cm high. The 
experiments were made with air flow at ambient conditions and with a 
centerline speed of 11.2 m/sec. 

In Fig. 3.4 the values of U (the streamwise mean velocity) obtained 
from a pitot probe and from the triaxial probe are compared. The pitot 
probe was modified to create the same stem blockage effect in the tunnel 
as did the triaxial probe, as seen in Fig. 3.5. In this test the probe 
was set to zero roll (a = 0°) and the axis was aligned with the flow 
(w = 0°). The readings of the pitot probe were corrected for shear dis- 
placement effect and for turbulence level. The maximum difference between 
the pitot probe and the triple wire probe occurs near the wall — about 
2.4%. The difference diminishes rapidly as the distance from the wall 
increases. The difference near the wall may be due to the finite size of 
the triaxial probe, as was explained earlier, interacting with the veloc- 
ity gradient. 

The effect of rotation of the probe around its axis was investigated. 
When the probe axis is aligned with the flow direction (m = 0°), rota- 
tion around its streamwise axis should not affect the result if the vel- 
ocity is uniform, but may affect the result in a shear flow. To investi- 
gate this, the probe was aligned with the flow direction (w = 0°), and 
for each transverse position across the tunnel, the probe was rotated 
around its axis to the values of a = 0°, 90°, 180°, and 270°. This 
angle range covers the extreme positions for the wires and exposes dif- 
ferent configurations of wires to the shear at different angles. 

Figure 3.6 shows the three mean velocity components (U, V, and W) 
as a function of the distance from the wall for several values of the 
roll angle a. Roll around the probe axis does not affect the U values. 
The effect on V and W is small, but not negligible. In this figure 
V and W should be zero, but due to the probe size, some deviation from 
zero is observed within the shear region. The most meaningful comparison 
for error in V and W is to compare them to the U at the same loca- 
tion. The largest deviations occur at the point near the wall for a - 0, 
V/U = 1%, and for a = 270°, W/U = 4.5%. These deviations from zero 
become smaller as the distance from the wall increases. In the zero gra- 
dient region at the centerline there is no deviation. 
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Figure 3.7 shows the turbulent kinetic energy and shear stress dis- 
tributions as a function of the distance from the wall for several values 
of the roll angle. The shear stress (-u ? v') measurements are compared 
with the linear shear stress distribution obtained from the pressure gra- 
dient along the channel (dP/dx) . As is seen, all the experimental data 
lies inside the ± 10% error band, but at angles a = 90° and a = 270° 
the deviations are much smaller; therefore, one would like to measure 
-u'v’ at these angles. The measured turbulent kinetic energy is not 
much affected by the roll. The largest difference between results occurs 
at the high shear region near the wall; it is about 3.5%. In the zero 
shear region there is no effect of roll. 

2 

Figu re 3.8 shows the diagonal Reynolds stress components (u’ , 

2 2 

v T , w 1 ) as a function of the distance from the wall, normalized with 

2 

centerline velocity. The streamwise normal stress (u* ) does not seem 

to be affected much by roll around the probe axis, even in the high shear 

2 

regions near the wall. On the same figure the u' distribution obtained 
with a conventional single horizontal wire is also given, and ins agree- 
ment with the triaxi al wi re data is good. The other normal Reynolds 
2 2 

stresses (v T and w’ ) are affected by the roll angle, especially in 
the high shear region, but the data collapse on each other quickly as the 
shear decreases. In the zero shear region on the centerline there are no 

deviations. On e important point to observe is that, at a certain a 

2 2 
value, if v ' reads high compared to the value at a = 0°, then w T 

2 

reads low, or vice versa, while u' does net change much with a. 

2 

This combination leads to q values which are quite insensitive to the 

2 

changes in a, a fortuitous result for the measurement of q — the 
main interest of the general research. 

The data discussed above are enough to qualify this system for mea- 
surements when the probe axis is aligned with the flow direction. But 
one of the most important objectives of this research was to find a method 
which would work in a flow of unknown direction. To investigate this, 
the probe axis was tilted against the approaching flow direction (w) , 
again in the 2-D channel. Some rotations around the probe axis (a) 
were also tested to see the combined effect of both a and w. The re- 
sults are discussed in the following paragraphs. 
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Figure 3.9 shows the three mean velocity components as a function of 
distance from the wall for several values of the angle between the flow 
and the probe axis (co ) . There the value of a was held constant, be- 
cause it was seen above that the mean velocities were not much affected 
by the roll engle. Up to oi = 20° the data for U collapse on top of 
each other and the deviation for w = 30° is not very large. The larg- 
est deviation at this angle is about 3.5% in the high shear region near 
the wall, and about 2% in the zero shear region. Deviation is calculated 
as the difference between two extremes, not from the pitot probe data. 
This result means that mean velocity can be measured with good accuracy 
if the approaching flow direction is within ±30° of the probe axis; 
i.e., one does not have to know the flow direction better than within a 
cone of half apex angle 30° around the probe axis to measure the mean 
velocity with acceptable accuracy. (If V = 0.1 U this will give an 
angle of about ±6%.) As the angle between the flow direction and the 
probe axis increases, the errors in V and W also increase. Some of 
this may be the effect of the probe size, as was explained earlier, and 

some may be due to the probe blockage effect in the channel. 

In Fig. 3.10 the turbulent kinetic energy (TKE) and shear stress 
are plotted for several values of rol] angle (a) and pitch angle (m) . 

In the TKE plot the line at the center is faired through the data at 

w = 0°, a = 0° (this measurement should be the one closest to reality). 
The other lines denote the ± 10% and ± 15% error bands around the refer- 
ence. Again the data points converge rapidly as the distance from the 

wall increases. Deviations are much smaller in the zero shear region. 

2 

In the same figure also, the q distribution for ct = 90°, a) = 20° is 

shown to demonstrate the increase in deviation as oo increases. Another 

important point to observe from this figure is that for a = 90° and 

a) = 10° the data lie very close to the center profile. This shows that, 

depending on the quantity being measured, there are angle combinations a 

and co for which the measurement cone can be enlarged. For example, for 

a = 90° and w = 10° and 20°, it appears that even in the highest 
2 

shear region q can be measured within 12% inside a cone of 15° half 

2 

apex angle around the probe axis. The deviations in q are not like 
uncertainty scatters, but rather have a preferred direction. It may be 
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possible to devise a scheme to correct the data based on the first 
estimate of the flow direction, to improve the accuracy. In the shear 
stress part of Fig. 3.10, the straight line in the middle of the figure 
is the shear stress distribution obtained from the pressure gradient in 
the streamwise direction. The other straight lines are the boundaries 
for ± 10% and ± 15% error. Most of the data up to the angle co = 10° lie 
within 10% error band, except a few points near the wall for angles a = 
270° and a = 90°. Almost all the data, including <u = 15° and a = 0°, 
lie within the ±15% error band. In conclusion, it can be said that the 
shear stress -u'v T can be measured within 10% within a cone of half 
apex angle 10° and within 15% inside a cone of half apex angle 15°, ex- 
cept very near the wal] . 

2 

Figure 3.11 shows the normal Reynolds stress u 1 as a function 

of the dis tance from the wall for several pitch angles. It was shown be- 

2 

fore that u’ does not change much with roll angle a. For this reason, 
different a and m combinations are not included in the figure, and 
attention is given only to the w values. Up to m = 15° the measure- 
ments collapse on top of each other except at the point nearest the wall, 
where the data for to = 15° deviate from others about 5%; but these data 
converge rapidly to the others with increasing distance from the wall. 

The data for to = 20° deviate about 15% from the data for smaller go's 

near the wall, but again conver ge fast with increasing distance from the 

2 

wall. The conclusion is that u T can be measured with good accuracy 

(maximum error being 5% near the wall) up to the angle oj = 15°, i.e., 

within a cone of half apex angle 15°. The deviation in the zero shear 

region is about 1% or 2% up to w = 30°. 

2 2 

Figure 3.12 shows the normal Reynolds stresses v’ and w f for 

2 

different roll and pitch angles. The data for v ? collapse on top of 
each other for w = 0° and w = 5°. For (w = 10°, ct = 0°) and ('jj = 
10 °, a = 180°), the data generally lie in the 10% error band, except the 
next-to-wall point for (w = 10°, a = 0°), for which the deviation is 
about 15%. Generally, it can be said that within a cone of half apex 
angle 5°, v ? can be measured with very good accuracy (less than 2% 
deviation) . The measurement cone angle can be enlarged to 15° in a zero 
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within 5% 


shear region. In the high shear region one can measure v' 
for the case of o> = 10°, a = 0°, and within 10-15% for oo = 10°, a = 

180 °. 

2 

Considering now the w T data, in the zero shear region the data 
converge very well, and good accuracies can be obtained up to w — 15°. 

Up to (0=5° the data lie within 10% of each other in the high shear 
region. For the larger values of w, certain combinations of a and aj 
are necessary to get the data within 10-15% error bands in a high shear 
region. For example, all of the u> = 10° data is out of the 15% error 
band for all a values, but the data for a; = 20° and a = 90° fall 
within the 15% error band. 

The qualification data discussed above show that the size of the 

measurement cone changes , depen ding on the quantity measured. Some quan- 

2 2 

tities, such as v’ and w* , are affected by rotation around the 
probe axis, especially in high-velocity gradient regions. The reason for 
these changes may be the large probe size. 

An important point shown by the qualification data is that the major 
mean velocity component U can be measured quite accurately within a 
cone of 30° half apex angle around the probe axis. In three-dimensional 
flows where the probe axis makes a large angle with the unknown flow di- 
rection, the errors in the measured turbulence quantities may be unaccep- 
table. In critical cases, one might use a two-step process, first finding 
the flow direction approximately (within 3-4°), and then aligning the 
probe with this direction. In this position all of the turbulence quan- 
tities could be measured quite accurately. 

3.7 Further Work on Triaxial Wire System 

As has been seen, the size of the probe creates some problems in 
high shear flows, by causing artificial velocity components to be indi- 
cated which depend on the arrangement of the wires in the velocity gra- 
dient. Work is in progress currently to predict the magnitude of the 
artificial velocity components indicated, as a function of the orienta- 
tion of the probe in flow and the value of the gradient of mean velocity. 
It is hoped that this work will lead to a method of correcting the turbu- 
lence quantities. 
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Secondly, to eliminate the probe size effects, a new, smaller probe 
design is being considered for possible future extensions of this work. 


3 . 8 Conclusions and Capabilitie s 

A hot-wire method and system have been developed to make measure- 
ments in 3-D turbulent flows with high fluctuations and partially unknown 
flow direction. The results are summarized below, The error figures are 
for high-velocity gradient (about 1600 1/sec) regions; for zero or low 
shear regions the errors are smaller and consequently the measurement 
cones are larger. With this system: 

• The mean velocity can be measured with good accuracy within a cone 
of 30° half apex angle around the probe axis. 

• The turbulence kinetic energy can be measured with 10-15% accuracy 
within a cone of half apex angle 10-12°. 

• The shear stress (-u f v f ) can be measured within ± 10% accuracy 

in a cone of half apex angle of 10° and ±15% accuracy within a cone 
of 15° half apex angle. 

• The streamwise normal Reynolds stress component can be measured 
with a maximum error of 5% within a cone of half apex angle 15°. 

The errors in the quantities mentioned above are not random but 
rather have a preferred direction. Some of these effects can be predic- 
ted from the known effects of probe size and position in the velocity 
field. Work in this area is continuing. 

The triple wire method is very fast, with the data-acquisition time 
being governed by the time required to establish a stationary value of 
the process being studied. Thus, large amounts of data can be taken. 

With the addition of RMS and correlating circuits, quantities such 
as TKE and mean velocity could be displayed directly with analog instru- 
ments. This would eliminate the need for complicated data reduction — 
usually the most tedious part of an experiment. 
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-D TURBULENT FLOW ANALYZER 
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Flow diagram of the complete system 









U (m/sec) 


Fig. 3.4. Comparison of triaxial probe and pitot probe 
measurements of U in the 2-D channel. 
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Fig. 3.5. Photograph of the triaxial probe and modified pitot probe. 




o 



q 2 (m 2 /sec 2 ) 


Fig. 3.7. Effect of roll angle about the probe axis on 
shear stress and turbulent kinetic energy. 
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turbulent kinetic 





2 

Fig. 3.11. Effect of pitch angle on the u’ 
Reynolds stress. 
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Chapter 4 


PREDICTION OF RECOVERY REGION HYDRODYNAMICS 
WITH A ONE-EQUATION MODEL OF TURBULENCE 


One of the main objectives of this research was to see whether or 
not the recovery region hydrodynamics could be predicted using a one- 
equat.ion model of turbulence in a two-dimensional boundary layer compu- 
tation program, STAN5 (Crawford & Kays, 1975). 

An algebraic equation for the mixing length was needed for the 
model. The actual mixing length profiles were first obtained from the 
experimental profiles (u’v ? and U) and purely empirical curve fits were 
used to model the behavior. The piecewise model was used in STAN5, as 
a check, and the predictions of the data were very successful. Each 
region in the mixing-length profile was then interpreted physically and 
alternative, plausible, equations identified. This procedure showed that 
the recovery region hydrodynamics can be satisfactorily explained by 
postulating a two-dimensional boundary layer growing inside the thicker 
initial boundary layer. The following sections present the details of 
these steps. 


4 . 1 Equations to be Solved 

The following equations must be solved to obtain mean velocity and 
turbulent kinetic energy (TKE) profiles. For a two-dimensional, turbu- 
lent boundary layer under isothermal conditions, at ambient temperatures 
and low speeds, at the constant free stream velocity and without any ex- 
ternal body forces: 

T T r) V 

Continuity: — b -r— = 0 (4.1) 

dx 9y 


x -Moment um : 



+ pV 


3U 

3y 


_ 3 _ 

3y 


(•* 


9U 

9y 


- p u ' v T j 


(4.2) 
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TKE: 


pu + pv 



V 


3y 


(J) (4.3) 


Boundary conditions: 
Momentum eqn. 


y 


) D 

It 


0 

0 


lim U 


y-Kn 


u 

oo 


(4.4) 


The equation for TKE is not solved all the way to the wall in STAN 5 , 

but only to y + = 2A + where A + is a measure of sublayer thickness. It 

is assumed that the flow is in local equilibrium (experiments confirm this 

+ + 

assumption) below y = 2A and Prandtl's mixing length can be used in 
this region. The TKE at this point is calculated such that at y = 2A 
the eddy viscosity obtained from the mixing length model is equal to the 
one obtained from a one-equation model of turbulence. This condition 
gives the following boundary conditions on TKE: 


2 

= 

2 




+ 

at y 


2A 


+ 


2 

lim 


2 

*eo 

2 


( 4 . 5 ) 


The following terms must be modeled in equations (4.2) and (4.3) in 
order to obtain a soluble set of equations: u ? v' , V , J . The term 
for u’v' will be modeled efter Boussinesq (1877) with eddy viscosity 
model: 


- u v 


= e. 


3U 
M 3y 


(4.6) 


The term 


'M 


will be modeled after Prandtl (1945) and Kolmogorov (1942): 


e 


M 



( 4 . 7 ) 


The term "V" dissipation of TKE will be modeled as given in Launder & 
Spalding (1972 ) : 
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V 


B K 

q 


(4.8) 



and diffusion of TKE "J M is also to be modeled as in Launder & Spalding 
(1972) . 


J 


(v + e q ) 


3(q 2 /2) 

9y 


(4.9) 


where Sc = £ /£ . The production of TKE, "P", 
q M q r * 


P = - P 


^g - »(f)*i?(g> 


can be expressed as 
- 2 


(4.10) 


To complete the model, the following constants or functions need to 
be specified: Z. A , B , Sc . A is the production constant and can 

q q q q 2 

be obtained from the value of the correlation coefficient - u’v’/q 
near the wall. It is generally about 0.22 (Wolfstein, 1969). is 

called the dissipation constant and can be evaluated under the condition 
that the production of TKE is equal to its dissipation in the region near 
the wall. This condition gives the following relation: 


B 

q 



(4.11) 


k = 0.41 


(4.12) 


From this equation, B = 0.377 is obtained. 

q 

"Sc^", the Schmidt number of TKE, is expressed as is shown in the 
following figure. It has a value of 1.75 near the wall and 0.5 near the 



0.0 0.5 1.0 

y/S .90 
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free stream. Launder & Spalding (1972) suggest that for film-cooling 
applications Sc^ should follow a linear distribution from 1.75 near the 
wall down to 0.5 near free stream, but during the present predictions, 
the distribution shown in the preceding figure was found to work better. 
The extended region of Sc^ = 0.5 near free stream played an important 
role in modeling the correct diffusion of TKE near the edge of the momen- 
tum boundary layer. The value near the wall did not seem to be as impor- 
tant because of the dominant role of production and dissipation in this 
region . 

The last quantity to be modeled to complete the set is the mixing 
length, Predictions and data show that all other constants men- 

tioned above have their usual values (i.e., the values used for 2-D flat- 
plate boundary predictions) . The success of predictions depended very 
strongly on the correct modeling of as will be shown and discussed 

in the following sections. 

4 . 2 Mixing-Length Model 

A mixing-length model was developed specifically for the recovery 
region. This same model can later be applied to the full-coverage re- 
gion, however, because the boundary layer should behave almost like a 
recovery region between the injection rows. 

The following requirements should be met by the mixing-length model 
for it to have at least some limited universality: 

• It should be possible to relate the deviations from the 2-D mixing 
length to physical events taking place in the flow field. 

• The dynamics of the model should allow it to relax back to a 2-D 
flat-plate mixing length. 

The general approach taken in modeling the mixing length was ex- 
plained at the. beginning of this chapter. More specifically, the model- 
ing process will be explained in Sections 4. 2. 1-4. 2. 9 in the following 
order: First the flow structure in the recovery region will be explained 

briefly (Fig. 4.1). Then each region in the mixing length profiles will 
be discussed and empirical equations will be given with supporting physi- 
cal arguments. Figures 5.27, 5.28, 4.1, and 4.2 are used for these 

discussions, and they will not be mentioned separately in the discussions 
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of each region. Figures 5.27 and 5.28 show the mixing-length distribu- 
tions obtained from experiments by 



(4.13) 


Three streamwise positions were used: the start of the recovery region 
(x = 188 cm) and two downstream recovery region stations (x = 214 cm, 
256 cm, or 27 hole diameters and 67 hole diameters downstream of the last 
row of injection. The mixing length at the start of recovery region was 
obtained from the spanwise-averaged profiles at this location. Fig. 4.2 
shows the modeled mixing length superposed on Fig. 4.1 to show the rela- 
tion between the regions of mixing length and the flow structure of the 
recovery region. The new information is on Z vs. y coordinates. 


4.2.1 Flow structure in the recovery reg ion 

Figure 4.1 shows a sketch of the flow structure in the recovery re- 
gion. The basic structure is an internal two-dimensional (2-D) boundary 
layer growing inside the thick boundary layer. The thickness of the 
inner 2-D boundary layer is 6' and the thickness of the overall bound- 
ary layer thickened with the injection process in the full-coverage re- 
gion is 6. In the middle regions the two boundary layers blend (the 
cross-hatched area) with the help of cumulative jet spread. The region 
next to the wall is not cross-hatched, however, because the processes in 
this region are completely controlled by the wall. The 2-D internal 
boundary layer has an initial thickness at the start of the recovery re- 
gion which depends on the upstream conditions (e.g., the blowing ratio). 
This internal boundary layer growth has been observed by several other 
experimenters whenever there is a sudden change of the surface conditions 
For example, Antonia 4 Luxton (1972) observed such an internal boundary 
layer in their experiments on the response of a turbulent boundary layer 
to a step change in surface roughness. 

When the mixing-length profiles are examined in Figs. 5.27 and 5.28 
five distinguishable mixing length regions can be identified, numbered in 
Fig. 4.2. In the same figure, three distinguishable flow regions are 
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shown, numbered I, II, and III. Region I is the near-wall region of the 
2-D internal boundary layer, II is the blend region dominated by cumula- 
tive jet effects, and region III is the outer region, dominated by the 
outer region of the thick boundary layer. 

4.2.2 Region 1 of the mixing length 

A region very close to the wall, where £ = <y, is termed region 1. 
Experiments show that this region extends up to y/<5 - 0.14 for M = 0.4, 
but only up to y/S - 0.055 for M = 0.9. 

This is the inner region of the 2-D internal boundary layer where 
the length scale is based on the distance from the wall "y". For pro- 
gramming convenience, the mixing length and the distance from the wall 
in all the regions were normalized on the total boundary layer of thick- 
ness 6. In this region, however, the wall effects are dominant and de- 
termine the heat transfer regardless of the blowing ratio. The effect 
of blowing is mainly to change the initial thickness of the internal 2-D 
boundary layer. The proposed model for this region is: 

(!) - -(!)» » < (!) - Kid «•>*> 

where D is the damping function (Van Driest) and (y/^^ i s t * ie depar- 
ture point from "tcy" line. This point corresponds to y/6’ = A/k for 
the inner 2-D boundary layer of thickness 6 1 . 

The following empirical equation is given for (y/S)^: 


(£) = A + 

C.M 

/ > 

| - C 0 

V 6 / k 

1 

\ 6 1 

' 2 

d 



— 


where = 0.0045 and ^2 ~ * It seen that as M gets large the 

initial value of (y/&)^ - A/k - C-j^M gets smaller, indicating that a 
higher blowing ratio destroys the near wall layer, causing a smaller ini- 
tial thickness for the 2-D internal boundary layer. This fits the physi- 
cal situation very well. C ^ is the number of boundary layer thicknesses 

at which (y/S). reaches the A/k point for the outer boundary layer, 
cl * 

In fact, then, the recovery is completed. (Antonia also supports C ? 
being 37.) controls the rate of recovery of (y/6)^ point. 

* Personal communication. 
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4.2.3 Region 2 of the mixing length 

Experiments show that in this region the mixing length is constant 
and low, compared to the outer region of a 2-D flat-plate boundary layer 
of the same thickness. This region is barely visible for M = 0.4 but 
is obvious for M = 0.9. It extends up to about (y/6) - 0.2 for all 
stations and corresponds to the outer region of the inner 2-D boundary 
layer, where £ = A 5’. Here the eddy size does not depend on the dis- 
tance from the wall, nor does it depend on the total thickness of the 
outer boundary layer. 

The distance of this region from the outer edge means that the outer 
length scale does not affect it. This region does not extend all the way 
up to y = 6 f for the reasons stated below. 

The following empirical equation is proposed: 



This region does- not extend all the way up to 5 ’ . As soon as the edge 
of the inner layer of the outer boundary layer (y / 6 = A/k) is encoun- 
tered, the length scale begins to be affected by the outer length scale 
of the outer boundary layer. The branch point is given as A/k. 

4.2.4 Region 3 of the mixing length 

Experiments show that in this region the mixing length is still below 
a 2-D flat-plate value. This is barely visible for M = 0.4 but is ob- 
vious for M = 0.9. In this region, the mixing length rises from the 
value in region 2 to 0.085. The region starts around y/6 = A/k - 0.2 
and extends up to y/6 - 0.3 for M = 0.9 and y/6 - 0.25 for M = 0.4 
for the first station. It moves away from the wall in the recovery re- 
gion. 

This is the first section of the blend region of inner and outer 
boundary layers (region II in Fig. 4.2). Here the length scale changes 
from the 2-D inner layer value to the value in the outer layer in a me- 
dium dominated by the jet flow regime. 

The following empirical equation is given for this region: 
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(f) = a (f) + b for 


< 


(I ) - (« )i 


(4.17) 


z 


The linear combination is all that is justified and fits the experiments. 
Here, (y/S)^ is the intersection point of the new mixing length line 
and the (Jl/6) = A line. The (y/<5)^ point is the effective center- 
line of rising jets and moves out in the recovery region as the jets 
spread outwards. This will later be shown quantitatively in section 
4.4.2. The following empirical equation is given for this point: 


/ y \ A ^ ^ 

/ x T \ 

Uk ~ k + c 3 m 

\~r) + l 4 _ 


(4.18) 


It has an initial value of A/k + C 3 C^M and an advance rate of C^M(x’/^) 
in the downstream direction. The initial value increases with M, indi- 
cating a deeper penetration for high blowing, which fits the physics. 

The values of the constants are: 


= 0.0275 , C, - 4.0 

3 4 


(4.19) 


and the equations for coefficients a and b in equation (4.6) are: 


a = 


A - k ( y/ 6 ) 


(y/6)i - A 


(4.20) 


b = 


- — + K( y /6 ) d 

(y/« ± - 


(4.21) 


4.2.5 Region 4 of the mixing length 

Experiments show that a region of augmented mixing length is located 
between 0.24 <_ (y/<5) < 0.37 for M = 0.4 and between 0.35 <_ (y/6) < 
0.55 for M = 0.9 at the start of the recovery region. This region is 
wider and higher for M = 0.9 than for M = 0.4. The augmented region 
vanishes in the recovery regions. 

The augmentation in the mixing length occurs due to high shear be- 
tween the cumulative jet spread and the outer layer fluid. The maximum 
in the augmentation occurs somewhere between the outer edge of the jet 
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spread and the effective jet centerline. In the case of higher blowing, 
this region is moved outwards because of the greater jet penetration. 

The augmented region also moves outwards in the downstream direction in 
the recovery region, as jets rise; but at the same time it vanishes, as 
there is a decrease in shear between the cumulative jet spread and the 
outer layer. This is discussed further in section 5.9. 

This region was modeled as a parabola, the simplest curve which fit 
the data: 


and 





(A. 22) 



(A. 23) 


Here, c, d, and e are constant coefficients. Their values depend on 
end points and will be given after the discussion of the end points. 

W is the width of the augmented region at the (£/<5) = A level. Based 
on the experimental evidence, it was made a function of blowing ratio M 
Since no better formulation was possible, it was expressed as a linear 
function. 


W = C 7 M (A. 24) 

In Figure A. 2, (Z/6) is the maximum value of the augmentation. Its 

° max 

decay in the downstream direction was expressed as an exponential func- 
tion. 

|) .(A) exp(-^p-) (4.25) 

max max , l ' b ' 

(Z/8) . is the initial value of augmentation. It is given from the 

max ,i 

experimental evidence, as follows: 

(!) = C 5 M (4.26) 

max, l 
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The coefficients for the parabola are given as: 


c = 


4 ( 4 / 6 ) 


W 


max 


8(Jl/6) 


d = 


max / v \ , W 

Is). + 2 

W 1 — 1 


e - - 


4 (4/6) 


W 


max 


m.Rn 


+ W 


The values of the constants are 


C 5 = 0.496 , C 6 = 0.435 , C y = 0.333 


( 4 . 27 ) 


( 4 . 28 ) 


(4.29) 


(4.30) 


4.2.6 Region 5 of the mixing length 

Experiments show that in this region the length scale is constant 
and it is (&/< 5) = A = 0.085. The region begins after the augmented re- 
gion and continues up to the. free stream for both blowing ratios. 

This region is the outer layer of the outer boundary layer, where 
the mixing length scales on the total thickness 6 . In this region, the 
following empirical equation is given: 



x for (f). +w < (f) 


(4.31) 


4.2.7 Summary of the model equations 


Region 1 (4.14): 


Region 2 (4.16): 


Region 3 (4.17) 
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Region 4 (4.23) (f) = * + (!) - (f). < (!) ± (!). + 

a i i 

Region 5 (4.31): (§ ) = * > (!). +w < (!) 


Unknowns in the above equations are: 
Region 1 (4.15): 


(!) - 7 +c i m [< x ’ /6) - c 2_ 


D = 1.0- exp (-y + / A + ) 


Region 3 (4.20): 


(4.21): 


(4.18): 


Region 4 (4.22) 


(4.27) 


(4.28) 


(4.29) 


(4.25) 
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(!). 
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A - k ( y/ 6) d 

(y/6) i - A/k 


- — + < (y/ 5) d (y/<5). 
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8(5,/ 6 ) 
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d = - 
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4(5,/ 6 ) 


max 


W 
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(!).[(!). + f 

i *— x — 1 


(!) = (!) . exp (- ^cf 1 ) 

max max,i ' 6 


(4.26) 


(!) 


max , x 


= C r M 
5 


(4.32) 
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( 4 . 24 ): 


W 


C y M 


Constants: 

= 0.0045 (controls the rate of recovery of (y/S)^)- 

= 37.0 (number of boundary layer thicknesses at which (y/6)^ 

approaches (A/k) or the 2-D internal layer coalesces with 
outer layer) . 

= 0.0275 (controls the rate at which the (v/<5)^ point advan- 

ces, i.e., the rising of the effective jet centerline). 

= 4.0 (controls the initial value of the (y / 6 ) ^ point at the 

start of the recovery region) . 

C,. = 0.496 (controls the magnitude of the initial augmentation). 

= 0.435 (decay rate of augmentation). 

C y = 0.333 (controls the width of the augmented region). 

Figure 4.3 shows the experimental and modeled mixing lengths for 
M = 0.4 and M = 0.9 at the three different stations in the recovery 
region. 


4.2.8 Comments on the number of constan ts 

Except for the usual universal values of sc = 0.41 and A = 0.085, 
seven constants were used to model the mixing length. The number of con- 
stants looks very large at first glance, but when examined more closely 
one can see that this number is reasonable. Three different flow struc- 
tures are being modeled: inner and outer boundary layers and blend re- 

gion in between, as well as their dynamics in the streamwise direction. 

During the prediction process it was observed that the most important 
constants in the mixing length model were C^, and C^. The con- 
stants used to specify the augmented region (C^, , and C y ) did not 

prove significant. The augmented region did not have much effect on the 
predictions because it lies in the outer region of the boundary layer. 

It also decays rapidly. During computer experiments with W and 
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(£/< 5) , changing the values of C.- and C-. did not affect the results 

at all. It is possible that the augmented region can be completely elimi- 
nated for the recovery region predictions, thus reducing the number of 
constants to four. The constants were obtained from the empirical data 
rather than from computer experiments. Kacker & Whitelaw (1970) used five 
constants to model the mixing length for prediction of wall jet and wall- 
wake flows, which is similar to, but not more complicated than, the re- 
covery region of the present full-coverage film-cooled surface. In the 
present model the augmented region is kept for convenience and for easy 
adaptation of this mixing length model to the full-coverage region where 
the peak of augmentation moves closer to the wall and is important for 
the prediction of heat transfer. 

In the recovery region predictions, especially for M = 0.9, the 
most important region of the mixing length turned out to be the reduced 
mixing length region near the wall. Without correct modeling of this 
region predictions always failed. The first four constants are the im- 
portant ones in the recovery region predictions. 

4.2.9 General remarks about the mixing length 

It was observed that only regions 2, 3, and 4 show deviations from 
a usual 2-D mixing length model. Region 2 proved especially interesting, 
because it had been thought that the jet mainstream interaction would 
increase the mixing length above the 2-D value (Choe et al. , 1975; Craw- 
ford et al., 1976). Choe et al . (1975) and Crawford et al. (1976) cal- 
culated the mixing length for low blowing and in the full-coverage region. 
They observed that the augmented region (Region 4) moved closer to the 
wall and did not observe Regions 2 and 3, which are more easily observable 
for high blowing. In reality. Regions 2 and 3 are very important in pre- 
dictions, as is explained in Section 4.3.1. 

4.2.10 Physical explanation of recovery to 2-D state in the mixing- 
length model 

The relaxation to the 2-D state takes place in the following manner. 
As the cumulative effects of jets move downstream and diffuse out, they 
lose their strength and the augmentation dies out ( 3 or 46), but the 
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effective centerline of the jets (the (y/S)^ point) continues to rise. 

At the same time the inner 2-D boundary layer thickness 6* grows faster 
than 6 (6 almost stops growing after the injection stops), and fi- 

nally, when they merge (in approximately 406), the recovery process is 
complete. 

4. 3 Predictions 

The mixing length model, along with the TKE equation, was used as a 
one-equation model of turbulence in the 2-D boundary layer program STAN5 
(Crawford & Kays, 1975). Starting with the spanwise-averaged mean veloc- 
ity and TKE profiles at the beginning of the recovery region (x = 188 cm) , 
very successful predictions of TKE and the mean velocity profiles for 
M = 0.4 and M = 0.9 were obtained at two downstream stations in the 
recovery region (x = 214 cm, 256 cm) . 

Figure 4.4 shows a comparison of the predicted mean velocity and TKE 
profiles with the experimental data for M = 0.4 on the two stations in 
the recovery region, starting with the spanwise-averaged initial profiles. 
The prediction of mean velocity on plate 11 (x = 214 cm) is somewhat 
high in the middle region of the boundary layer, but the difference les- 
sens by the 27th plate (x = 256 cm). For TKE profiles the predictions 
are somewhat low. These figures show that the model predicts a slightly 
faster recovery to the 2-D state than the physical process, but the dif- 
ference is not that great. The reason for lower TKE predictions might 
be the mixing length. In Fig. 4.3 the mixing-length model for M = 0.4 
is slightly lower compared to measurements in the region 0.15 < (y/6) < 
0.3. This means that the model results in a smaller production and a 
higher dissipation than the reality, yielding lower TKE values compared 
to the experiment. 

Figure 4.5 shows a comparison of the mean velocity and TKE profiles 
with the experimental data for M = 0.9 on the two stations in the re- 
covery region, starting with the spanwise-averaged initial profiles. The 
suggested model predicts the M = 0.4 case well. The predictions for 
the M = 0.9 case are excellent. It is important to remember that the 
constants do not change with the blowing ratio. Both cases were predic- 
ted with the same set of constants for the mixing length and the other 
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constants in the one-equation model are the usual ones employed for 2-D 
boundary layer predictions. The author believes that the model will suc- 
cessfully predict blowing ratios up to M = 1.0, Further extrapolation 
of the model might be dangerous due to the changing hydrodynamic charac- 
ter of the flows. 

4.3.1 The effect of reduced mixing on predictions 

As previously mentioned, it is important to consider the effect of 
the reduced mixing near the wall in making predictions. This is demon- 
strated by Fig. 4.6. This figure shows the comparisons of the experimen- 
tal data with the predictions of TKE (in the recovery region on the 27th 
plate for M = 0.9). The predictions were made with two different mixing 
length models. One of them is the normal 2-D flat-plate mixing length, 
and the other is the model developed in this study, which has a smaller 
value near the wall relative to the former. The new model predicts per- 
fectly, whereas the usual 2-D mixing length predicts much higher TKEs 
almost up to the first half of the boundary layer. The reason lies in 
the relationship between the mixing length and the production and dissi- 
pation of TKE. Higher mixing lengths increase the production and reduce 
the dissipation, giving rise to higher TKEs. 

The reduced mixing region did not have a significant effect on near- 
velocity profiles; however, it i_s_ important in heat transfer behavior 
because of its influence on the turbulence level. 

4 .4 Refining the Model 

The empirical curve fits described in section 4.2.7 do in fact match 
the boundary layer behavior. The physical arguments are plausible and 
suggest that the empirical relations could be replaced (if desired) by 
more conventional forms. In this section some possible alternative equa- 
tions are examined — without changing the physical arguments — both for 
the purpose of testing the physical argument for reasonableness of the 
magnitudes and to take advantage of established forms where possible. 
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4.4.1 The outer edge of the inner boundary layer 
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Using a conventional equation to describe the growth of the internal 
boundary layer growth (Schlichting, 1968) yields: 
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Therefore, 
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Consequently, if one knows the initial value of 6 and (y/<S) at 

d 

x f = 0, then the initial 6 f can be calculated and x* determined 

o 

from the formula. After x’ is found, the rest of 6' and (y/5) , 

o d 

can be obtained from the boundary layer growth formula. 

In Eqn. (4.33) the value of U can be used for the value of U 
since U at is never far from U^. For example, for M = 0.9 

the velocity profile is so flat that near 6’^ 

U ^ U 


For M = 0 . 4 the initial thickness of this boundary layer is large (be- 
cause there is less disturbance), and again around y 't 6 ' , U ^ U . 

One might question the use of Eqn. (4.33) to predict the growth of the 
inner layer, since it neglects the effect of the turbulence of the outer 
regions. It seems justified, however, because around 6 ’ the turbulence 
levels are quite small for both cases studied here (3%-5%) . 
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By using U - = 16.7 m/sec, p =1.2 kg/m , y = 0.9 x 10 

kg/msec, the virtual origin of the internal 2-D boundary layer can be 
calculated for both cases from the initial value of (y / 6 ) ^ - The follow- 
ing values were obtained: 


For 

M = 0.4 , 

x T = 1.29 m 

0 

For 

M = 0.9 , 

x ’ = 0.41 m 

o 


(x^)^_q ^ is greater than 9* is normal because for low 

blowing, the undisturbed region near the wall is larger. After calcula- 
tion of virtual origins, the following comparisons were obtained between 
the piecewise model and the internal boundary layer growth formula. 


M = 0.4 


(x’/S) 

6’(m) 

0 

^ 5 

-v 12 

Piecewise Model 

0.028 

0.036 

0.045 

Boundary Layer 
Formula 

0.028 

0.033 

0.040 


M = 0.9 


<S’(m) 

0 

^ 5 

^ 12 

Piecewise Model 

0.011 

0.018 

0.028 

Boundary Layer 
Formula 

. 

0.011 

0.017 

0.025 


As is seen from the comparison, the present linear model and the boundary 
layer growth results compare quite well. In fact, when one looks at the 
data, the boundary layer formula seems to be closer. 

The following changes can be suggested in the present model. Instead 
of calculating from the piecewise linear formula (Eqn. (4.15)), 

(y/6) d can b e calculated, which involves only one empirical con- 

stant. Then, using 
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Other 6 f can be calculated as follows: 


and 


/ P U (x'+x')\" 1/5 
’ = ( x * + x q) 0.37 ( jp- 2-) 

(!) - (!) (x) 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


Here, even though the equation for boundary layer growth is still empiri- 
cal, it is more universal than the constant supplied in the formula for 
(y/<5)^. So the (y/6)^ formula can be reduced to one empirical constant 
supplied by the present experiments and a more universally accepted em- 
pirical equation. 


4.4.2 Comparison of the development of the (y/5)^ point with jet 
spreading theory 

It is argued that the dynamics of the (y/6)^ point corresponds to 
the rising of the effective jet centerline. In this section, the output 
of the empirical equation for (y/<5)^ (Eqn. (4.18)) is compared with 
the results of an. equation describing the rising of the centerline of a 
jet in crossflow (Abramovich, 1960) : 



Here x” is the virtual origin for the effective jet action, D is the 
diameter of a single jet, and a f is the injection angle. The solution 
of this equation depends on four parameters: x", D. M and a’. The last 

three are supplied by the problem physics, and the only unknown is the 
virtual origin, which can be calculated from the empirical input, i.e.. 
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from the initial value of (y/6)^, assuming that effective jet spread 
occurs according to the equation (4.41), 

The results of the virtual origin calculation are: 


For 

M = 0.4 , 

x" = 
o 

0.13 m 

For 

M = 0.9 , 

x" = 
o 

0.055 m 

Then, with these results 

and from Eqn. 

(4.41) 

, the spreading of the jets 


can be calculated and compared with the present model values of (y/fi)^- 
The results: 


M « 0.4 


^^\(x’/6) 

(y/<s) i 

0 

'v- 5 

^12 

Piecewise Model 

0.244 

0.3 

0.37 

Jet Spread For- 
mula 

0.244 

0.3 

0.333 


M = 0.9 


(y/ <5 ) ± 

0 

v 5 

^ 12 

. 

Piecewise Model 

0.3 

0.424 

0.59 

Jet Spread For- 
mula 

0.3 

0.55 

0.67 


As is seen again, the comparison is not bad. This supports the argument 
that (y/<5)^ point moves out with distance similar to the jet centerline 
The following changes can be made in the model to replace the con- 
stant with the more accepted empirical equation; or with this physi- 

cally supporting agreement, the old model for (y/S)^ can be left as it 
is. The change is as follows. Calculate 
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Then 
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(4.44) 


4.4.3 Other changes suggested 

To include the effect of P/D (pitch- to-diameter ratio, the blowing 
ratio M in the formulas can be replaced by 

F = 

or by 

F' = 

which reflects the relative effects of the jet and free stream momentum. 
It is thought that it will be better to replace them by F T ; by inclu- 
sion of an empirical jet spread equation (4.10), the variation in the 
injection angle (a’) can be represented in the model. 
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2-D MIXING LENGTH 
MODELED MIXING LENGTH 



(y/8) d (X/k) (y/8)j y/8 99 1.0 


Fig. 4.2. Mixing-length model and flow structure in the 
recovery region. 





Fig. 4.3. Comparisons of the piecewise mixing-length 
model and the data in the recovery region. 
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PREDICTION 

DATA 



Predictions of streamwise mean velocity component (U) and turbulent 
2 

kinetic energy (q ) for M = 0.4. 







Fig. 4.6. Effect of mixing-length on predic- 
tions of turbulent kinetic energy. 
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Chapter 5 

MEASUREMENTS AND DISCUSSIONS OF THE DATA 

Measurements were made to supply the necessary empirical input for 
development of turbulence models for the recovery region following a 
full-coverage, film-cooled region. The final objective of the turbulence 
modeling was to obtain correct predictions of the heat transfer coeffi- 
cient. This chapter will discuss the experimental conditions, measure- 
ment locations, and the data. Due to the large volume of the data, only 
the parts which directly concern the prediction process will be discussed 
in detail. 

5 . 1 Experimental Conditions 

The experiments were made under isothermal conditions, with air tem- 
peratures around 22-25°C. Since the heat transfer experiments were made 
with small temperature differences (16°), it was believed that the hydro- 
dynamics of the isothermal flow field would be the same as that of non- 
isothermal heat transfer experiments. This question was investigated by 
Harnett, Birkebak & Eckert (1961). The pressure levels were about atmos- 
pheric 760 mm of mercury) , and the experiments were made at a uniform 
free stream velocity of approximately 16 m/sec. The velocity of the free 
stream was kept uniform along the test section by adjusting the top wall, 
(as explained in Crawford et al . , 1976), so that the static pressure 
change along the tunnel was less than ±0.8% of freestream dynamic head 
(± 0.005 inches of water). 

Experiments were made at two different blowing ratios: M = 0.4 and 

M = 0.9 (M = U /U ) . The M = 0.4 blowing ratio was chosen because, 
jet °° 

during heat transfer experiments, minimum Stanton numbers on the test 
section were obtained here, as demonstrated by Crawford et al . (1976). 

The second blowing ratio, M = 0.9, was chosen because a definite change 
in the heat transfer behavior was seen at this high blowing ratio. 
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5 . 2 Measurement Locations 

The schematic of the test section is given in Fig. 2.4. The jets 
are injected into the main stream at a 30° angle with the horizontal 
plane, and the injection is in line with the main flow direction. Due 
to the interaction of the jets and the main flow, a 3-D high-fluctuation 
flow field is created, at least locally, around the secondary air jets. 
All the measurements of the velocity field were made with the 3-D Turbu- 
lent Flow Analyzer system introduced in Chapter 3. 

The measurement locations can be seen in Fig. 5.1. Five streamwise 
positions were selected for measurements to observe the development of 
the flow field on the full-coverage surface and its relaxation in the re- 
covery region. The first is located at the center of the 5th plate (x = 
148 cm) after three rows of blowing, and the second is at the center of 
the 9th plate (x = 168 cm.) after seven rows of blowing in the full- 
coverage region. The third one is at the center of the first plate (x = 
188 cm) in the recovery region after 11 rows of blowing. The last two 
are at the centers of the 11th (x = 214 cm) and 27th plates (x = 256 cm) 
and are, respectively, 27 and 67 hole diameters downstream from the last 
row of injection holes. In order to see spanwise variations in the flow, 
several measurements were made at each streamwise position. At the first 
and third stations, five spanwise locations were used, situated at z/P = 
+ 05, + 0.3, 0.0, -0.3, and -0.5. Since the flow field is periodic, it 
repeats itself after z/P = + 0.5 and - 0.5. At the second streamwise 
location three spanwise locations were selected at z/P = 4- 0.5, 0.0, 
and - 0.5. In the recovery region for M = 0.9, at each streamwise 
location three spanwise locations were used. After the experiments it 
was seen that the profiles in the recovery region station were the same 
in all three spanwise locations, as illustrated by Figures 5.4 and 5.3. 

So for M = 0.4 these spanwise locations were eliminated in the recovery 
region. The profiles were used to obtain the necessary spanwise-averaged 
input to the 2-D boundary layer program STAN5. Measurements were made at 
each location of the three mean velocities and the six Reynolds stresses. 
Each quantity was averaged over 22 seconds to include the effects of low- 
frequency f luctuations . 
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5 . 3 Mean Velocity Profiles 

In this section attention will be given to the streamwise component 
of the mean velocity U, rather than discussing the three-dimensionality 
of the flow around the holes. Two sets of profiles for different blowing 
ratios will be discussed and compared with each other. 

5.3.1 Mean velocity profiles for M = Q.4 

Figure 5.2 shows the mean velocity profiles at each spanwise loca- 
tion, at each of the five streamwise stations, for the blowing ratio 
M = 0.4. The location of each of the traversing stations is marked on 
each figure with respect to the injection locations. 

The first observation on this figure is the remarkable symmetry of 
the profiles from positions which are symmetric with respect to the cen- 
terline. This proves that there is periodicity in the spanwise direction 
and that the coalescing of the upstream jets is regular. Periodicity is 
very helpful in spanwise averaging. 

In the full-coverage section, two distinct regions appear in each of 
the profiles. In the outer region all spanwise profiles are the same; 
however, in the near-wall region, differences are observed in the span- 
wise profiles because of injection. The same observation was made by 
LeBrocq, Launder & Priddin (1971) for a 45° injection. The spanwise pro- 
files al so show some common behavior. The central profiles jusc down- 
stream of a jet show the greatest velocity defect, because they are 
affected directly by injection. The profiles at the outer edges show 
less defect than the central profiles but more defect than the intermedi- 
ate positions (between two holes). This is expected, because the outer 
profiles are affected by the injection one row upstream and are still 
recovering from the effect of that injection, while the profiles between 
two holes have no injection in line with them and are affected only by 
the lateral spreading of jets. This is clearly seen in Fig. 5.2, where, 
after three rows of blowing, the intermediate profiles show almost no 
defect, thus indicating that the effect of lateral spreading has not 
reached them with its full force yet. After 11 rows of blowing, however, 
the momentum defect in the intermediate profiles can easily be observed. 


73 


The size of the near-wall momentum defect region increases down- 
stream. The velocity defect region after 11 rows of blowing is almost 
twice as large as after three rows of blowing. This indicates the cumu- 
lative effect of upstream jets which keep rising as they spread down- 
stream. The jets penetrate only about 2 cm, immediately after injection, 
as can be seen from the figure, the point where there is an abrupt change 
in the slope of the. profiles. The immediate penetration distance is about 
two jet diameters for this blowing ratio. This means that the mass shed 
into the boundary layer by each injection row penetrates to an almost 
constant distance for a certain blowing ratio at the injection location 
before jets are knocked down by the main stream. The immediate penetra- 
tion distance increases slightly in a downstream direction as the bound- 
ary layer loses momentum due to the effect of previous injections. This 
constant penetration distance was also observed by Le Brocq, Launder & 
Priddin (1971), Crawford et al. (1976), Colladay 6 Russel (1975), and 
Ramsey & Goldstein (1971). 

There is no indication of a reverse flow and separation of the jet 
from the surface in any of the full-coverage region profiles. This might 
be due to the fact that the profiles were measured at locations where 
jets have already reattached themselves to the surface. The measurement 
locations were 2.5 hole diameters downstream from the edge of the closest 
hole. The separation of jet from the surface was observed by Bergeles, 
Gosman & Launder (1975) for normal blowing at a blowing ratio of M = 0.24, 
and extended only to 1.5 diameters downstream of the hole. It is gener- 
ally agreed in the literature (Le Brocq, Launder & Priddin, 1971; Colla- 
day & Russell, 1972) that slant- injection jets stay attached to the sur- 
face at much higher blowing ratios compared to a normal injection. 

Launder & York (1973) state that for a 45° in-line injection with P/D = 

8, jets attach to the surface at 3-4 hole diameters downstream of the 
last injection for M - 0.6. These observations help one arrive at the 
conclusion that the jets for the 30° slant injection for M = 0.4 re- 
attach to the surface in a very short distance downstream or perhaps 
don't separate at. all. 

The points which are closest to the wall in all the mean velocity 
profiles have almost the same value of U/U^ for the same spanwise 
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locations at all streamwise stations. For the central profiles the value 
is U/U^ = 0.35-0.4. This number can be obtained by multiplying the 
value of the blowing ratio (th^/U^ = 0.4) with the cosine of the in- 
jection angle (30°). This is a definite indication of the average jet 
velocity in the near-wall region. It was observed by Bergeles, Gosman & 
Launder (1975) that the velocity in the exit plane of jets is not uniform 
at all; however, the present work shows that even though the velocity is 
not uniform, the effect is governed by the average jet velocity given by 

M x U . 

00 

th 

In Fig. 5.2, curves representing (1/7) power profiles are shown 
at the start of the recovery region (after the last row of injection) and 
at the two measurement stations in the recovery region as reference lines 
to make the cross-comparisons easier. It was expected that, finally, in 
the recovery region, the boundary layer would return to a normal 2-D 
flat-plate boundary layer. A (1/7) ^ power profile was chosen to rep- 
resent the flat-plate mean velocity profile, although any curve which is 
close to a 2-D turbulent boundary layer profile could have been chosen 
as a reference line. A spanwise-averaged profile is also shown (for the 

profiles at the start of the recovery region only) . When compared to 

th 

the (1/7) power profiles, it can be seen that the momentum defect de- 
creases as the boundary layer recovers from the blowing. However, there 
is a considerable momentum deficit even at the last recovery region sta- 
tion due to the large initial deficit. The reason for the long recovery 
distance is the "momentum- sink" action of the wall. The diffusion of 
momentum from the free stream is the only way for these profiles to re- 
cover to 2-D flat-plate profiles, and diffusion is a slow process. The 
long recovery distance indicates that the effects of jets near the wall 
prevail over long distances for low blowing. 

Velocity gradients at the wall are smaller than for a 2-D flat-plate 
profile; this will give lower friction factors. There is also no indica- 
tion of boundary layer separation except locally at the injection points. 
Thus, the aerodynamic behavior of the 30° slant-hole in-line injection 
with M = 0.4 is quite good. Since this blowing also gives the lowest 
heat transfer coefficients (Crawford et al., 1976), it is a very important 
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blowing ratio in full-coverage film-cooling applications. Le, Brocq, 
Launder & Priddin (1971) agree with this conclusion. 

Even though the jets initially penetrate to a constant distance (two 

ttl 

hole diameters), comparisons with the (1/7) power profiles show that 
their cumulative effect on the momentum deficit can be observed at much 
larger distances from the surface due to the rising of the effective jet 
centerlines in the downstream direction. The jets continue to rise even 
at large distances from the last row of holes. This observation was also 
made by Keffer and Baines (1963) . 

Figure 5.3 shows the spanwise- averaged mean velocity for M = 0.4 
in (y/6) and U/U^ coordinates. The profile at x = 168 cm is not 
to be trusted quantitatively near the wall. Only three spanwise profiles 
were taken at this location, and the spanwise averaging of three profiles 
does not represent a physically correct average. This point is further 
explained in section 5.6. It is plotted here, however, for comparisons 
in the outer layers where the spanwise profiles are not very different 
from each other, and the resulting average is correct. The profile after 
three rows of blowing (x = 148 cm) has the largest momentum defect, but 
the difference from the others is not great. This means that the region 
near the wall is governed by the latest injection and stays almost at the 
same level of U/U^. The effects of previous rows can be seen in the 
outer layers. For example, in the outer regions, the profile at x = 

148 cm is very close to the profile at the last recovery station (x = 

256 cm). This means that the jets cannot strongly affect the outer 
layers at this distance. The outer region profiles are very similar after 
7 and 11 rows of blowing (x = 168, 188 cm, respectively) and in the re- 
covery region profile at 27 hole diameters downstream cf the last row of 
injection (x = 214 cm) . The result of these observations can be summa- 
rized as follows: after four or five rows of blowing it appears that the 
boundary layer reaches a state of equilibrium in a large part of the 
boundary layer (almost above y/6 >_ 0.1) for low blowing. This may be 
very helpful in the modeling of the injection process in the full- 
coverage region. 
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5.3.2 Mean velocity profiles for M = 0.9 

Figure 5.4 shows the spanwise mean velocity profiles at five stream- 
wise stations for the blowing ratio of M = 0.9. Spanwise symmetry can 
be obs'erved, as in M = 0.4. This indicates uniform lateral jet spread. 
Nina & Whitelaw (1971) cautioned against non-uniform coalescing of jets in 
full-coverage film-cooling applications; such a phenomenon is not observed 
here. 

Again, two distinct regions can be observed in the profiles, an outer 
and an inner one. In the outer region all the spanwise profiles are the 
same at a certain streamwise station. In the inner region, however, there 
are differences in the profiles, depending on the spanwise location. The 
size of the inner region is very large compared to the M = 0.4 case be- 
cause of the higher rate of blowing and deeper penetration of the jets. 

The jets penetrate up to three hole diameters at M = 0.9, and this dis- 
tance is almost the same for all the stations. 

Another difference between the cases is the excess momentum near the 
wall in M = 0.9 because of the high blowing ratio. This can be observed 
by comparing the spanwise-averaged profile to the (1/7) ^ power reference 
curve at the third streamwise location. 

In the central profiles just downstream of a jet, a local maximum in 
velocity occurs around y = 1 cm with U/U = 0.85-0.9 (this was not 
observed for M = 0.4). This point is the centerline of the new injected 
jet, about one hole diameter above the surface. This can be observed at 
the same location in every station of the full-coverage region. This 
maximum in the streamwise mean velocity was also observed by Le Brocq, 
Launder & Priddin (1971) for M = 0.5 and P/D = 8 for 45° injection. 
This is another indication of the fact that the depth of immediate pene- 
tration is constant. 

It is important to observe that at the point closest to the wall in 
the central profiles, there is a sharp decrease in the velocity. This 
indicates that the jet has separated from the surface and has not reat- 
tached yet at the measurement location (about 2.5 diameters downstream 
of the hole). However, the same phenomenon is not observed in the side 
profiles, which are 7.5 hole diameters away from the last injection lo- 
cation. Therefore, it can be concluded that, for the case of M - 0.9, 
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jets reattach to the surface somewhere between 2.5 and 7.5 hole diameters 

/ 

downstream of their injection location. This was also observed by ColTa- 
day & Russell (1975) and by Launder & York (1973) . The separation of the 
protective jet from the surface explains the poor heat transfer behavior 
(high heat transfer coefficients) of the M = 0.9 case observed by Craw- 
ford et al. (1976) . 

In Fig. 5.4, (1/7) ^ power profiles are shown at the start of the 

recovery region and at the two recovery region stations as reference 
lines for comparison. A spanwise average is also shown for the profiles 
at the beginning of the recovery region. 

The spanwise-averaged mean velocity profile has a region of negative 
slope (as can be seen in section 5) which gives rise to a locally nega- 
tive shear stress. 

There is a rapid relaxation back to the 2-D state in the recovery 

region, as can be seen from the comparison of the profiles in this re- 
th 

gion with (1/7) power profiles. In fact, the profile at the last re- 

til 

covery region station is almost the same as the (1/7) profiles, except 
for a small momentum deficit in the outer layer. Relaxation back to the 
2-D state is much faster for M = 0.9 than for M = 0.4 because of the 
excess momentum near the wall. The wall acts as a momentum sink and fa- 
cilitates the relaxation process. In the outer layer, where there is 
some momentum deficit, the relaxation process is very slow. The momentum 
supplied by the free stream entrainment and by the transport is not large 
enough for rapid relaxation. 

Higher velocity gradients near the wall (when compared with the 

t - Vi 

(1/7) profile and the M = 0.4 profiles) lead to larger skin friction 
coefficients, thus increasing the aerodynamic drag. Increase in the 
skin friction coefficitn with blowing was observed by Kacker & Whitelaw 
(1970) for a 2-D wall jet. 

One observes generally the lack of velocity gradient all over the 
M = 0.9 mean profiles except at the wall. This yields low shear 
stresses and low TKEs (low turbulent mixing) in a large part of the 
boundary layer. Samuel & Joubert (1965) also observed that the closer 
the M is to unity, the lower the turbulent mixing is for a 2-D 
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injection. Seban & Back (1962) made the same observation in the case of 
a 2-D tangential injection. 

In Fig. 5.5, spanwise-averaged profiles for five streamwise loca- 
tions are plotted on (y/S) and U/U coordinates for M = 0.9. A very 
clear outer layer similarity is observed in these coordinates; the pro- 
files lie much closer to each other than in the same case for M = 0.4. 

The similarity is not observable below (y/6) = 0.25. The high rate of 
blowing causes non-equilibrium in the region near the wall. Another in- 
dication of this non-equilibrium can be seen in Fig. 5.4. The interme- 
diate profiles in each streamwise location get fuller downstream as the 
preceding jets expand laterally; they do not reach an equilibrium even 
after 11 rows of blowing. 

Figure 5.6 shows the spanwise-averaged mean velocity components in 
semi-logarithmic coordinates (U/U vs. Jin y) for M = 0.4 and M = 0.9. 
No logarithmic region can be seen in the profiles for either case, in the 
full-coverage region. The signs of the curvatures of the profiles in the 
full-coverage region are opposite for different blowing ratios; this is 
an indication of the momentum deficit in one case and of the momentum 
excess in the other. In the recovery region, however, there is a defi- 
nite indication of a logarithmic region near the wall. The two profiles 
at x = 214 cm and x = 256 cm exhibit linear regions for both blowing 
ratios; this proves that as soon as the blowing stops, wall effects are 
dominant and the "law of the wall" seems to apply. 

5.4 Turbulent Kinetic Energy Profiles 

Turbulent kinetic energ} T will be discussed, like the mean velocity 
profiles, in the form of two sets — first for M = 0.4, second for 
M = 0.9. Comparisons between the two cases will be deferred until the 
section discussing the M = 0.9 data. 

5.4.1 Tur b ulent kinetic energy profiles for M - 0.4 

Figure 5.7 shows the TKE profiles at several spanwise locations at 
each of five streamwise stations for M = 0.4. Spanwise symmetry can be 
observed in these profiles, indicating an orderly behavior of the turbu- 
lent flow field. 
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The position of the 6 ^ thickness is marked on the profiles in 

the full-coverage region. As seen, the TKE boundary layer extends up to 
y - 1.3 6 This was also observed by several other workers — for 

example, Klebanoff (1955). 

In the full-coverage region, there is a strong convergence above 
y = 2 cm in the TKE profiles at all spanwise locations, but large dif- 
ferences exist between the profiles below y = 2 cm. This relates to 
the 2 cm penetration distance and indicates that the region near the 
wall is locally affected by the jets. 

Near the wall, the highest TKE is observed in the central profiles; 
the second-highest one is observed in the side profiles. The intermedi- 
ate profiles have the lowest TKE, as expected (no jets are in line with 
them) . There is a dip in the center profiles around y = 1 cm. This 
can be explained as follows. There is a high velocity gradient near the 
wall which causes high shear, resulting in higher TKEs. However, there 
exists a very low velocity gradient region (see Fig. 5.2) just above it, 
because of injection. In this region the local TKE production is de- 
creased but local dissipation is high, leading to low TKE values. The 
low velocity gradient region ends at about y = 2 cm (penetration dis- 
tance). Above this level, higher velocity gradients lead to higher TKE. 
As a result, a double maximum is observed in the central TKE profiles. 

The side TKE profiles have only one maximum. At some distance from the 
injection location, the TKE near the wall dissipates and the maximum near 
the wall vanishes. TKE diffuses from the second maximum towards the wall. 
The result is a smooth profile with one maximum. 

In Fig. 5.7, 2-D flat-plate TKEs are also shown both at the start of 
the recovery region and at the two recovery region stations as references 
for comparisons. The 2-D flat-plate TKE profiles were obtained with the 
STAN5 computer program (Crawford & Kays, 1975) for the same free stream 
turbulence level as the experimental profiles. A spanwise-averaged TKE 
profile is shown for the start of the recovery region. Comparisons with 
the 2-D profile show that there is excess TKE throughout the region, but 
most of the excess energy lies in the middle section of the boundary 
layer. There is more violent turbulent mixing in the outer layers rela- 
tive to layers next to the wall. Even though the jets do not penetrate 
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very far into the boundary layer, the TKE diffuses towards the free 
stream, causing large TKE values far from the wall. 

In the recovery region, there is still high TKE compared to the 2-D 
value. The reversion of the turbulence quantities to the 2-D state is 
very slow. This interesting process is shown in the following sketch. 



Again, there are two distinct regions in the relaxation process, an outer 
region and an inner region. The relaxation to the 2-D state is in the 
expected direction in the outer region and is caused by the diffusion of 
TKE. However, an opposite "relaxation" is seen in the inner layer. 

This means that the rate of diffusion of the TKE toward the wall 
exceeds the rate of dissipation near the wall. After the profile 
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becomes flat in the last recovery region station (67 diameters downstream), 
the diffusion will stop and the dissipation will dominate, causing the re- 
version of the profile to the 2-D state. 

The high TKE (turbulent mixing) levels in the recovery region must 
be the reason for the slowly rising Stanton numbers here for M = 0.4. 
Crawford et al. (1976) observed that Stanton numbers for M = 0.4 de- 
crease in the full-coverage region and rise slowly in the recovery re- 
gion. This process can be explained by the TKE levels, as follows. The 
TKE levels are higher than the 2-D boundary layer values. This higher 
turbulence level should cause high Stanton numbers; however, at the same 
time, an energy sink is created in the full-coverage region by the cold 
injection. The combination of these two opposing phenomena results in 
the low and decreasing Stanton numbers in this region. There is no in- 
jection in the recovery region (i.e., the energy sink is lost), but the 
turbulence level (turbulent mixing) is still high. The mixing process, 
unopposed, increases the Stanton numbers and they rise gradually in the 
recovery region. 

I n Fig . 5.8, spauwise-averaged TKE profiles are plotted in (y/5) 
and Vq^/U coordinates for each streamwise station. There is no defi- 
nite similarity in any region, even though all the profiles seem very 
similar (again, for the x = 168 cm profile, the points near the wall 
are to be neglected, as discussed in Sec. 5.6). There is an indication 
of invariance in the TKE profiles in the full-coverage region after four 
or five fows of blowing (over (y/6) = 0.0). This was also the case with 
the streamwise mean velocity profiles. 

5.4.2 Turbulent kinetic energy profiles for M = 0.9 

Figure 5.9 shows the TKE profiles for several spanwise locations at 
each of five streamwise stations for M = 0.9. 

Spanwise symmetry can be observed in the profiles within the experi- 
mental uncertainty limits. 

The behavior of the TKE is qualitatively similar to the M = 0.4 
case — contrary to the mean velocity profiles which exhibit different 
behavior. Once more, the highest TKE is seen in the central profiles, 
the next highest in the side profiles, and the lowest in the intermediate 
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profiles. However, quantitatively, there are big differences between the 
TKE profiles for the M = 0.4 case and the M = 0.9 case. 

The double maximum observed in the M = 0.4 TKE profiles is not 

found here. The highest TKE occurs next to the wall rather than in the 
outer region, owing to the high velocity gradient there. Above y = 

2.0 cm there is a sharp drop in TKE. Also, the TKE levels for M = 0.9 

are much lower than in the case of M = 0.4 within a large region of 

the boundary layer outboard of y/6 = 0.15 (compare the spanwise-averaged 
profile at the start of the recovery region with the 2-D TKE profile) . 

This is clear evidence of a decrease in turbulence level as the blowing 
ratio approaches 1. The TKE levels are lower because of the lack of 
velocity gradients for the M = 0.9 case, as shown in Fig. 5.11. This 
figure shows the central profiles of mean velocity and TKE at the start 
of the recovery region for M = 0.4 and M = 0.9. The mean velocity 
gradient is smaller and the TKE is lower for the M = 0.9 case than for 
the M = 0.4 case everywhere above (y/6) ~ 0.2. This behavior may be 
due primarily to the dependence of the TKE production on the mean velocity 
gradient — the larger the mean velocity gradient, the larger the TICE pro- 
duction and the TKE itself. 

The lower TKE levels for the higher blowing case seem, at first 
glance, to be inconsistent with the heat transfer data obtained by Craw- 
ford et al . (1976), but when examined closely it can be seen that there 

may be an explanation: the governing process for the heat transfer may 
be different. Crawford et al. (1976) observed that the Stanton number 
for the M = 9.9 case is below the 2-D flat-plate value and above the 
value for the M = 0.4 case. The Stanton number values decreased in the 
downstream direction in the full-coverage region and also in the recovery 
region. If turbulent mixing were the governing process in the heat trans- 
fer, one would expect the Stanton numbers for the M = 0.9 case to be 
lower than for the M =0.4 case, for two reasons: first, more cold 
fluid is injected into the boundary layer for high blowing and, second, 
there is less turbulent mixing. Instead, higher Stanton numbers were ob- 
served. The reason could be that, for high blowing, the jets separate 
from the wall and penetrate deeply into the boundary layer, entraining the 
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hot free stream fluid towards the wall in the middle lanes. This process 
is schematically explained in the following figure. 



This figure was also used by Le Brocq, Launder & Priddin (1971) to de- 
scribe the in-line blowing process. This same process seems to occur in 
a staggered array, at high blowing. Thie hypothesis is also supported 
by the V component of the mean velocity, as can be seen in Figs. 5.16 
and 5.17. 

In addition to the entrainment process, the cold fluid (energy sink) 
is also deposited in the outer layer due to the separation of the jets 
from the surface. Thus, two opposing mechanisms result in the decreas- 
ing Stanton numbers in the full-coverage region. (They are, nevertheless, 
higher than the low-blowing values.) 

In summary, the two pairs of opposing processes which control the 
heat transfer are: (i) turbulent mixing versus energy sink near the wall 

for low blowing, and (ii) outer-layer fluid entrainment versus energy 
sink in the outer layer for high blowing. 

In the recovery region, the Stanton numbers continue to decrease 
because there is rapid decay of entrainment of the outer-layer 
fluid by the jets and the energy sink is no longer being injected 
at the wall. In the absence of these mechanisms, the process continues 
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in its natural direction until the turbulent mixing is strong enough to 
reverse the trend. Since the three-dimensionality is important in heat 
transfer for high blowing, any 2-D computer program will have difficulty 
in predicting the heat transfer in the full-coverage region for these 
cases . 

In Fig. 5.9, the recovery region TKE profiles are compared with 2-D 
flat— plate TKE profiles. The behavior of the TKE is very interesting. 

At the start of the recovery region, the spanwise-averaged TKE is higher 
than the 2-D value; however, by 27 hole diameters into the recovery re- 
gion, TKE values are lower than the 2-D value below y - 3 cm. In fact, 
the deviation from the 2-D value increases in the downstream direction. 
Normally one would expect the difference to decrease. This can be ex- 
plained as follows. In Fig. 5.4, the spanwise-averaged mean velocity 
profile at the start of the recovery region has negative velocity gradi- 
ents close to the wall. There are also negative shear stresses in the 
same area as seen in Fig. 5.13. In fact, since the flow is three- 
dimensional, there are y locations where the mean velocity gradient is 
positive and the shear stress negative, and vice versa. This results in 
the apparent negative production of TKE. The high TKE near the wall is 
dissipated rapidly as soon as the injection stops and thus cannot feed 
the low TKE regions by diffusion. As a result of the ''negative produc- 
tion", the dissipation, and the decreased diffusion, the TKE falls below 
the 2-D flat-plate value. Even though the negative velocity gradient 
region diminishes rapidly after the last row of injection, its effect on 
TKE continues because the turbulence profiles respond slowly to changes. 
The production first becomes zero and then positive. A considerable dis- 
tance is required before the production can adjust itself to increase the 
TKE level. 

The relaxation of TKE in the recovery region occurs in the expected 
direction in the outer layer. This rate is again slow, as was observed 
in the M = 0.4 case. In the near-wall region, TKE responds rapidly. 

At the point nearest the wall , TKE drops from V?/». = 0.134 to 0.1 
quickly, after the last row of injection, and remains at that level 
throughout the recovery region. This is also true for M = 0.4, in the 
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recovery region (see Fig. 5.7) and for a 2-D flat-plate boundary layer. 
These observations illustrate the rapid recovery to the 2-D state in the 
near-wall layer (y/6 _< 0.1), which occurs as soon as the injection 

stops. 

Figure 5.10 shows the spanwise-averaged TKE profiles for M = 0.9 
at the five streamwise stations (discard the x - 168 cm, as discussed 
in Section 5.6). There is a fair similarity above (y/6) = 0.8. This is 
due to the undisturbed TKE region here. Below (y/6) = 0.8, there is 
similarity between the full-coverage region profiles and the recovery 
region profiles. This may mean that a state of invariance for TKE is 
reached in the full-coverage region. This phenomenon may simplify pre- 
dictions of the full-coverage region and might lead to better insight 
into the injection process. The similarity observed in the recovery re- 
gion may be due to the slow rate of relaxation to the 2-D state. 

5.5 Streamwise Shear Stress Profiles 

Figures 5.12 and 5.13 show the streamwise shear stress profiles at 
the start of the recovery region (just after the last row of injections) 
for M = 0.4 and M = 0.9, respectively. The 2-D flat-plate shear 
stress profiles (Klebanoff, 1955) are also shown in these figures as 
reference lines for comparison. In both cases there is some symmetry ob- 
served for the profiles which are in symmetrical positions with respect 
to the centerline. Compared to the 2-D shear stress, excesses are ob- 
served for the M = 0.4 case, above y/6 - 0.2, for all of the span- 
wise profiles. In general, for the M = 0.9 case, however, defects can 
be seen when compared to the reference profiles below y/6 ~ 0.6. Above 
y/6 - 0.6, there are some excesses, but the magnitudes of these are less 
than for the M = 0.4 case owing to the higher shear action of the ef- 
fective jet spread for the low blowing. Thus, the shear stress levels 
are much higher for M = 0.4 than for M = 0.9, as was expected from 
the mean velocity gradients. This is another indication of high turbu- 
lence mixing for the low blowing. The extremely high shear stress next 
to the wall in the central profile for the M = 0.9 case must be due to 
the separated jet, which creates a high-shear and vortex region there. 
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Some negative shear stresses can be seen for M = 0.9 in the cen- 
tral (No. 3) and side profiles (Nos. 1 and 5) due to the separation of 
jets. These negative shear stress regions occur close to and at the 
points of negative streamwise mean velocity gradients. It appears that 
-u’v’ is more or less following the direction of oU/3y. Shear stress 
regions close to zero can be seen in the intermediate profiles (Nos. 2 
and 4) next to the wall. These regions correspond to the 3U/3y = 0 
region of the mean velocity profiles. 

Another indication of the jet separation from the surface for M = 
0.9 is the highly disturbed region of shear stress in the central pro- 
file extending up to (y/6) - 0.5 (about y - 3 cm or three hole diam- 
eters). This is the penetration distance for the last row of injections. 
In fact, the same disturbed region can also be observed in the side pro- 
files because of the jets injected one row before. However, the profiles 
in the middle lanes are very smooth and do not show any sign of distur- 
bance. They also indicate the low level of turbulence for the M = 0.9 
case. 

Figure 5.14 shows spanwise-averaged streamwise shear stress profiles 
at five streamwise positions for M = 0.4. The 2-D flat-plate shear 
stress profiles are also shown for comparison. The development of the 
shear stress in the full-coverage region and its relaxation in the recov- 
ery region can easily be observed. In the full-coverage region the pro- 
files grow fuller in the downstream direction as the injection effects 
diffuse outwards. In the recovery region the relaxation process is quite 
slow. The shear stress levels are much higher in comparison to the 2-D 
values at each station, indicating high turbulent mixing. 

Figure 5.15 shows spanwise-averaged streamwise shear stress profiles 
at the start of the recovery region, as well as at the two recovery re- 
gion stations for M = 0.9. As can be seen, the relaxation process in 
the recovery region is also slow for this case. When compared to the 
2-D values, the shear stresses are quite small, indicating low turbulent 
mixing below (y/6) - 0.65. In the outer layers (above (y/6) - 0.65), 
there is some excess shear, but it is not significant. The low shear 
level is consistent with the mean velocity gradients at the start of the 
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recovery region. The spanwise-averaged mean velocity profile is quite 
flat, resulting in small gradients and small shear (see Fig. 5.4). The 
mean velocity reverts back to the 2-D flat-plate value rapidly in the 
recovery region. However, because of the slow relaxation speed of the 
turbulence field, the shear stress remains low for long distances. 


5 . 6 On Spanwise Averaging 

Most of the data above was discussed on spanwise-averaged bases. 

This section will explain the hows and whys of spanwise averaging. 

In section 5.7 it is shown that the three dimensionality of the flow 
field was local and restricted to the vicinity of the jets. The flow 
field was expected to be periodic due to the fact that the jets are in a 
regular array. The spanwise uniformity and the symmetry of the flow 
field were shown in the discussions of the mean velocity, TKE , and the 
shear stress profiles. 

The values of the heat transfer coefficients sought here are not the 
local values , but rather the average behavior of the Stanton numbers in the 
streamwise direction. This requirement plus the properties from the pre- 
ceding paragraph allow for a simpler two-dimensional analysis instead of 
a 3-D one. 

The following method was used to spanwise average the flow proper- 
ties. A spline fit in the lateral (z) direction was made to the pro- 
files at each y level within a pitch. Each point of the spanwise aver- 
aged profiles was obtained from the following formulas. To be consistent 
with the physics, mean velocities and TKEs were mass averaged. 


G 

average 



U G dz 


U dz 


(5.1) 


where G stands for U or q , as the case may be. Shear stresses 
were area averaged. 
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As was seen in section 5.2, five profiles were taken at the first 
and last full-coverage stations and three profiles at the second station. 
Results of spanwise averaging showed that three profiles were not enough to 
obtain a good spanwise average, since this process gave more weight than was 
proper to the central profile. Therefore, for the profiles at x = 168 cm 
(second full-coverage region station) , little importance is attached to 
the results of spanwise averaging near the wall, when the spanwise pro- 
files differ significantly from each other due to the injection. How- 
ever, spanwise averages obtained from five profiles appear to represent 
the physics adequately. The spanwise-averaged profiles of mean velocity 
and TKE at the start of the recovery region bears a strong resemblance 
(at least from the standpoint of qualitative tendencies) to the profiles 
in the recovery region, where no spanwise averaging was done because all 
the spanwise profiles are the same. From this it can be concluded that 
the physical averaging process taking place in the flow must coincide 
with the spanwise averaging of five profiles here. 

5 . 7 On the Three-Dimensionality of the Flow Field 

As noted, the three-dimensionality of the flow field was limited to 
the close vicinity of jets; the flow had a strongly preferred direction. 
Elaboration follows in sections 5.7.1 and 5.7.2. Space and time did not 
permit plotting all the data (such as the components of the mean velocity 
and Reynolds stress tensor) at every measurement location. Some interest- 
ing representative data of this nature are given, however, as an aid to 
understanding the three-dimensionality and its implications for heat 
transfer . 

5.7.1 V and W components of the mean velocity 

Figure 5.16 shows the V component of the mean velocity after the 
last row of injection for several spanwise locations for the case of 
M = 0.4. Together with the V, the flow angle y is also shown. Only 
the profiles on one side of the centerline are plotted because of sym- 
metry. As expected, because of the inline jet, the highest V occurs in 
the centerline profile at the point next to the wall. The smallest 
V values occur in the intermediate profile, which has no jet in line 
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with it. It is significant that the direction of V depends on the pro- 
file position. At the central and side positions, V is in the positive 
direction (towards the free stream), whereas for the intermediate profile 
it is in the negative direction. The fluid is pushed upwards by jets at 
the position in line with them and is entrained towards the wall in be- 
tween. Maximum values of V occur near the wall, reaching approximately 
8% of the free stream velocity; whereas in most of the boundary layer 
they are only from 0-4% of the free stream velocity. The largest flow 
angle is about 10° and occurs next to the wall. The flow angles get 
smaller rapidly as the distance from the wall increases. 

Figure 5.17 shows the V component of mean velocity after the last 
row of injection at several spanwise locations, for M = 0.9. The obser- 
vations for M = 0.4 also apply here qualitatively. Quantitatively, the 
values of V are higher because of the higher rate of blowing — there 
is a larger drift towards the wall in the middle lanes. This is evidence 
of the strong role which entrainment must play in the heat transfer for 
high blowing ratios. 

The largest value of V is around 16% of the free stream. The 
largest flow angle is about 12°. 

Figure 5.18 shows the W component for M = 0.4 at the same stream- 
wise location as the V components. An additional spanwise profile on 
the other side of the centerline is also shown. W is biased, not symmet- 
ric. Except for a few points near the wall, the W values are insignifi- 
cantly small (1-3% of the U^) . The flow angles 6 are also plotted on 
the same figure. Generally, V and W values for the M = 0.4 case are 
small enough so that locating the spanwise-averaged flow T as a 2-D boundary 
layer is a valid approximation. 

Figure 5.19 shows the W component for M = 0.9. The profile loca- 
tions are the same as on Figure 5.18. There is no preferred direction of 
the W, as seen in the case of M = 0.4; indeed, the W values are 
smaller than the values for M = 0.4. This is because of the low span- 
wise gradients in the streamwise mean velocity as compared to the low 
blowing case. The lower W values support both the spanwise averaging 
and the weak three-dimensionality argument, except perhaps very close to 
the wall, where W values are about ± 2% of the TJ , 

7 OO 
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The large and highly disturbed values of V and W near the wall 
for M = 0.9 in the central profile are indications of vortices created 
by the separated jet. 

In both cases of blowing, the angle range of the flow fiel^ is en- 
tirely suitable for measurements, within a good accuracy, of most of the 
turbulence quantities with the triaxial wire system, except perhaps at 
the point next to the wall just downstream of a jet for the M = 0.9 
case. 

5.7.2 Reynolds stresses 

Figure 5.20 shows the Reynolds shear stresses for the first station 
in the recovery region (27 hole diameters downstream of the last row of 
injection) for M = 0.4 and 0.9. These profiles are given to show a 
sample of the data; the rest are found in Appendix a* These profiles are 
not discussed in detail because of their secondary importance in the pres- 
ent work. 

The (u’v*) component of shear stress for M = 0.4 is signifi- 
cantly higher than for the others. The (u'w') and (v'w') components 
are almost one order of magnitude smaller than (u’v'), except for a 
few points near the wall. This indicates a return to two-dimensionality 
for the stress field. 

All of the shear stress values for M = 0.9 are low in comparison 
to M = 0.4 as expected, with the exception of (u'v’) near the wall, 
which is high (still remembering its upstream history). 

Figure 5.21 shows the normal Reynolds stresses at the same station 
as in Fig. 5.19 for M = 0.4 and M = 0.9. Again, the general observa- 
tion is of the expected low turbulence level of M = 0.9 in the outer 
region. The two cases also differ in the near isotropy of the M = 0.9 
case. All the normal stress components are very close in magnitude, 
which is not the case for M = 0.4. This distribution of the normal 
stress may be attributed to the flat mean velocity profiles and the 

small spanwise gradients in the M = 0.9 case. 

2 2 

For M = 0.4 the u' profile looks very much like the u’ pro- 
file obtained by Wilson & Goldstein (1976) in a turbulent plane wall jet, 
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except that here the levels are lower. The reason is that the plane wall 
jet was injected to a still environment, thus causing higher shear rate. 

5 . 8 Correlation Coefficients 

Figure 5.22 shows the spanwise-averaged stress-energy ratio 

2 

(-u’v’/q ) for M = 0.4 at five streamwise stations. The 2-D flat-plate 
value (0.15) is also given in the same figure for comparison (taken from 
Hinze, 1975). The values near the wall are smaller than the 2-D values 
because of the disturbance caused by the injection there. The values 
reach the 2-D level above (y/6) = 0.5, where the jet effects are weak. 
Recovery to the 2-D values near the wall can be observed in the recovery 
region stations. 

The first profile after three rows of blowing does not show the large 
area of low stress-energy ratio shown in the other full-coverage region 
profiles. This must be because the effects of injection have not yet pen- 
etrated that far. 

Figure 5.23 shows the same profiles for M = 0.9. A much larger area 
of reduced and disturbed stress-energy ratio is observed, attributed to 
the deeper jet penetration; otherwise the same observations can be made 
for this case as in M = 0.4. 

Figures 5.24 and 5.25 show the spanwise-averaged 3-D stress-energy 
2 

ratio, ( | T | / p ) / q , for M = 0.4 and M - 0.9, respectively. Since 

there is some three-dimensionality in the flow field, it was thought that 

it would be interesting to see the effect of three-dimensionality. Fig. 

5.24 shows that the total stress-energy ratio lies closer to the 2-D value 
2 

than - u ' v T / q , except for a few points near the wall. This suggests 
that the value of 0.15 may not be restricted to 2-D boundary layers, but 
may be useful in weakly 3-D flow such as this one. The same observations 
can be made for M = 0.9, also, except that the effect of high blowing 
still shows itself in large disturbances. 

Figure 5.26 shows the correlation coefficient, ~u’v f / 
for both blowing ratios at the two recovery region stations. The flat 
plate value (0.45) of the correlation coefficient is also shown for 
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comparison (Schlichting, 1968) . Again, for low blowing, the values are 
much closer to the flat-plate value than for high blowing. The effect of 
the injection can be observed near the wall in the form of decreased cor- 
relation coefficients for both cases. Reversion to the 2-D flat-plate 
state can be seen in the downstream direction for both blowing ratios. 


5.9 Mixing Length 

The mixing length can be obtained using Prandtl's (1925) definition 


(4.13) 
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Figures 5.27 and 5.28 show the mixing length obtained from the shear 
stress and the mean velocity data using the above formula in the recovery 
region for M = 0.4 and 0.9. The same figures also show the mixing- 
length profile for a 2-D flat plate boundary layer (from Escudier, 1966). 
An augmented mixing-length region is observed at the station just after 
the last row of blowing (x = 188 cm) for both cases. This augmentation 
is a result of the injection process, which creates shear stresses as 
well as regions of low 9U/9y. Evidence in support of this hypothesis is 
found in the large region of augmentation for the M = 0.9 case, where 
the velocity profiles are quite flat. Having a high mixing length does 
not always imply a high rate of mixing; it may just be a result of defi- 
nition. For example, even though low shear stress and TKE levels are ob- 
served for M = 0.9 (Figs. 5.9 and 5.13), the height and width of the 
augmented region are much larger than for M = 0.4. This is an indica- 
tion of the inadequacy of the mixing-length models in regions with a zero 
mean velocity gradient. Fortunately, the augmented region moved out at 
the start of recovery region so that its effect on the near wall behavior 
is negligible. Furthermore, the peak in the mixing length dies out quite 
fast, and is not observable in the recovery region stations. 

Another interesting point is the existence of regions of low mixing 
length near the wall, as compared to the 2-D flat-plate case. This region, 
cannot be seen clearly for M = 0.4, but is quite clear for M — 0.9. 

The details of this region were discussed in Chapter 4, along with the 
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predictions. Except for these two regions, the mixing length behaves 
like a 2-D mixing length. 

Figure 5.29 shows the mixing-length distributions obtained from the 
spanwise-averaged profiles in the full-coverage region for both blowing 
ratios, as compared to a 2-D flat-plate boundary layer. Some very gen- 
eral remarks can be made on these trends, even though the whole picture 
is quite complex. In the outer layers, the mixing length is similar to 
that for a 2-D layer. In the inner layer, the profile for M = 0*9 is 
reduced below the reference. However, a different behavior can be ob- 
served for M = 0.4. It appears that the profiles at the first and sec- 
ond recovery region stations have higher values than the 2-D below 
(y/<5) -0.15, at which point their values drop below the reference be- 
tween (y/S) - 0.15-0.4. This leads one to the conclusion that for 
low blowing ratios, the peak in the mixing length occurs near the wall 
and moves outwards in the downstream direction in the full-coverage re- 
gion, whereas for high blowing, it occurs in the outer layers to start 
with, due to deeper jet penetration. The points which are outside the 
figure limits are shown on top of the figure, with arrows indicating 
their magnitude. 

Figure 5.30 shows the mixing length obtained from the TKE using the 
following definitions: 

- u'v' = e M |2 , e M = * Vq 2 /2 (5-3) 

It was interesting to observe its behavior, since a one-equation model of 
turbulence was used in predictions. Again, there is some scatter in the 
data; the same trends discussed before apply here. Therefore it was de- 
cided to use the mixing length obtained from Eqn. (4.13) — since it was 

already programmed in the computer program by Crawford & Kays (1975). 

Finally, as a point of interest, the shape factor H is shown in 
Fig. 5.31 for both blowing ratios. It was obtained from the spanwise- 
averaged profiles. A representative value of H for a 2-D boundary layer 
(by Hinze, 1975) is also shown for comparison. The shape factor lies very 
close to the 2-D value for both blowing ratios. This is an indication of 
spanwise-averaged aerodynamic behavior which is not violently damaged. 
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It is over the 2-D value for M - 0.4 and below it for M = 0.9. The 
shape factor converges towards the 2-D value downstream in the recovery 
region for both blowing ratios. 
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FULL-COVERAGE H— RECOVERY REGION 



96 


Fig. 5.1. Measurement positions 
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Profiles of the strearawise mean velocity component (U) 





Fig. 5.3. Spanwise-averaged profiles of the streainwise 
mean velocity component (TJ) for M = 0*4. 
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Fig. 5.7. Turbulent kinetic energy profiles for 
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Fig. 5.9. Turbulent kinetic energy profiles for 
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Streamwise turbulent shear stress component 
at the start of the recovery region for M ; 
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Streamwise turbulent shear stress component (-u T v T ) at the 
start of the recovery region for M = 0.9. 
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Fig. 5.14. Spanwise-averaged profiles of the streamwise turbulent shear stress 
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Fig. 5.15. Spanwise- averaged profiles of the streamwise turbulent 
shear stress component for M = 0.9. 




Fig. 5.16. V component of the mean velocity at 
the start of the recovery region for 

M = 0.4. 
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Fig. 5.17. V component of the mean velocity at 
the start of the recovery region for 
M = 0.9. 


Ill 


III 






-10 -5 0 


$ (degrees). 

Fig, 5.18. W component of the mean velocity 
at the start of the recovery re- 
gion for M = 0.4, 
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Fig. 5.19. W component of the mean velocity at 
the start of the recovery region for 
M = 0.9. 
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Fig. 5,20. Reynolds shear stresses at the first recovery 
region station. 
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Fig. 5.21. Reynolds normal stresses at the first recovery 
region station. 
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Fig. 5.22. Spanwise-averaged stress-energy ratio 
for M = 0. 4 . 
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Fig. 5.23. Spanwise-averaged stress-energy ratio 
for M = 0.9. 
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Fig. 5.26. Correlation coefficient for M = 0.4, 0.9 
at the two recovery region stations. 
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Fig. 5.30. Mixing-length obtained from the turbu- 
lent kinetic energy in the recovery region 





Chapter 6 

SUMMARY AND RECOMMENDATIONS 


6.1 Summary 

The work reported here is a part of continuing study of the heat 
transfer and hydrodynamics associated with full-coverage film cooling. 

In the present work, hydrodynamic measurements were made with a triaxial 
hot-wire in the full-coverage region and the recovery region of a 30° 
slant-angle injection pattern. These measurements were made in isothermal 
flows at ambient temperature and pressure (about 20°C and 760 mm Hg) for 
two blowing ratios: M = 0.9 and M = 0.4. Profiles of the three mean 

velocity components and the six Reynolds stresses were obtained at sev- 
eral spanwise positions at each of five locations down the test plate: two 
within the full-coverage region, two within the recovery region, and one 
at the dividing line between the two regions. 

One of the most important observations from these experiments is the 
peculiar dependence of TKE on the blowing ratio: high TKEs were observed 
for low blowing (M = 0.4), low TKEs for high blowing (M = 0.9). This 
phenomenon results from the low mean velocity gradients for high blowing 
(close to unity). The difference in the TKE levels has a significant ef- 
fect on heat transfer behavior. 

As a result, it was concluded that two different pairs of processes 
compete to govern the heat transfer, depending on the level of blowing: 
for low blowing, the energy sink near the wall competes with high turbu- 
lent mixing; for high blowing, the energy sink in the outer layer competes 
with the entrainment and convection of mainstream fluid toward the wall in 
the lanes between the jets. The present 2-D spanwise-averaged solution 
methods are not expected to predict the heat transfer well for high blow- 
ing ratios — except in a purely formal way -- because they cannot handle 
the entrainment process and by definition neglect the spanwise 3-D flows. 

The flow in the recovery region can be described in terms of a two- 
layer model: an outer boundary layer where the length scale scales on 
the total thickness of the layer, and an inner layer where the mixing 
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length is the usual icy and A6. The two layers blend into each other 
with the spreading of jets. Recovery to a 2-D boundary layer is completed 
when the inner and outer layers finally merge. 

A one-equation model of turbulence was used in the 2-D finite- 
difference boundary layer computer program (STAN5) to predict the mean 
velocity and TKE profiles in the recovery region. The one-equation model 
employs the TKE conservation equation with an algebraic relationship for 
the mixing length. Mixing-length values calculated from the data were 
input to the program using a piecewise continuous, heuristic fit consistent 
with the concept of the two quasi- independent layers observed in the re- 
covery region. This mixing length pattern, used with a set of otherwise 
normal constants (for 2-D boundary layer predictions), successfully pre- 
dicted all the spanwise-averaged features of the flow. This strongly sug- 
gests that the principal mechanisms were modeled adequately by the mixing- 
length formulation. It was also shown that the piecewise continuous heur- 
istic model can be replaced by a set of relations taken from the literature 
describing the spreading of jets and the growth of boundary layers, thus 
supporting the physical arguments behind the mixing-length model. 

Large amounts of data were required by the present program, and a 
new hot-wire scheme was developed for this work: a triaxial hot-wire with 
an analog device for real-time data reduction. The method is basically 
simple; i.e., with three orthogonal wires, the components of the instan- 
taneous velocity vector can be deduced without recourse to any time aver- 
aging and without invoking the low-fluctuation assumption. The present 
method can tolerate an unknown flow direction within the following limits. 
The mean velocity was measured within a cone of ± 30° half apex angle 
around the probe axis with a maximum error of 4%, even in a high-velocity 
gradient (about 1600 1/sec) . The turbulent kinetic energy was measured 
within a cone of 12° half apex angle with a maximum error of 12% up to 
the same velocity gradient point and within 5.5% in a region of zero vel- 
ocity gradient. The method is very fast and practical for taking large 
amounts of data in three-dimensional, high-fluctuation turbulent flows. 
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6 . 2 Recommenda tions 


• The present model for the length scale, developed for the recovery 
region, can probably be used also for the full-coverage region provided 
that the injection process is modeled correctly. 

• Either 3-D solution methods should be used or new 2-D solution pro- 
cedures should be devised to simulate the entrainment process for high 
blowing ratios, in order to better predict the heat transfer. 

• For the 3-D turbulent flow measurement scheme, a smaller probe 
should be manufactured to decrease the errors arising from the probe 
size and to permit a closer approach to the wall. 

RMS and correlator circuits can be added to the 3-D Turbulent Flow 
Analyzer which will increase the speed of data- taking, as well as allow- 
ing for reading important quantities such as mean velocity and TKE di- 
rectly from the dials. This will greatly increase the range of applica- 
bility. 
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APPENDIX A 


TABULATION OF EXPERIMENTAL DATA 


A.l. Individual Profiles 

This section contains the dimensional and non-dimensional profiles 
of three mean velocity components and six Reynolds stresses at each mea- 
surement location for blowing ratios M = 0.4 and 0.9. Profiles of tur- 
bulent kinetic energy, flow angles, mixing-length, and correlation coeffi- 
cients are also given. The number of diameters downstream at the top of 
appropriate pages represents the distance of the measurement location from 
the last row of injection in terms of jet diameter D. The nomenclature 
given in this section also applies in Section A. 2. 

Nomenclature 

BETA 3, angle, between U and W components of mean velocity in degrees 

CM cm . 


DELM 

<5, 

boundary layer thickness. 

DELM1 

6 r 

displacement thickness. 

DELM 2 

V 

momentum thickness. 

GAMA 

Y, 

angle between U and V components of mean velocity in degrees 

L 


mixing-length. 

M 

P. U. ^/Poo U m» blowing ratio, 
jet jet “ “ 

M/S 

m/ sec . 

M2/S2 

2, 
m / 

2 

sec . 

P 

P, 

hole spacing or pitch. 

PR 

P, 

pressure. 

Q 

V? 

- Vu- 2 + v' 2 +W - 2 

Q2 

2 

q 

2 2 2 

= u f + v’ + w' , turbulent kinetic energy. 

R 

R 

/ 2 2 2 

= YU + V + W , magnitude of mean velocity. 
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RF 

RF2 

TF 

U 

UDERV 

UF 

UF2 

UINF 

UP 

UP2 

UV 

UW 

V 
VP 
VP 2 
VW 


w 

WP 

WP2 

X 

Y 

Z 


U , value of R at freestream. 

CO 

—2 

R at freestream. 

T^, freestream temperature (°C). 

U, component of mean velocity in x-direction. 

9U/8y, gradient of U component of mean velocity. 

U , freestream velocity. 

2 

U , freestream velocity squared. 

U , freestream velocity. 

00 7 J 

u', component of fluctuating velocity in x-direction. 

2 

u’ , Reynolds normal stress in x-direction. 

u’v', longitudinal-normal velocity correlation, Reynolds shear 
stress . 

u'w', longitudinal- tangential velocity correlation, Reynolds 
shear stress. 

V, component of mean velocity in y-direction. 

v 1 , component of fluctuating velocity in y-direction. 

2 

v’ , Reynolds normal stress in y-direction. 

v'w', normal-tangential velocity correlation, Reynolds shear 
stress . 

W, component of mean velocity in z-direction. 

w’, component of fluctuating velocity in z-direction. 

2 

w f , Reynolds normal stress in z-direction. 

x, distance along test surface measured from nozzle exit. 

y, distance normal to test surface. 

z, distance in transverse direction on test surface measured 
from centerline. 
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M*0.4 BLOWING REGION AFTER 7 ROWS OF HOLES# X-160 CM AND Z/P**0.5 

RUN 12147601 TF* 21.00 C PR* 760.00 MM-HG DELM* 3.629 CM DELM1* 1.105 CM D6LM2* 0.629 CM 
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M= 0.4 BLOWING REGION AFTER 7 ROWS OF HDLESt X«I 68 CM AND Z/P« 0.0 

RON 12137606 TF* 20.50 C PR» 760.50 MM-HG OELM* 3.628 CN DELM 1 - 1.204 CM 0 ELM 2 - 0.614 CM 
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M»O.A START OF RECOVERY AFTER 11 ROWS OF HOLES, X-1B0 CM AMO Z/P* 0.0 

RUN 12137603 TF= 20.50 C PR=» 760.50 MM-HG DELM- A.A19 CM 0ELM1= 1.301 CM DELM2- 0.760 CM. 
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A . 2 . Spanwis e-Aver aged Profiles 

’ This section contains the spanwise-averaged profiles of streamwise 
mean velocity component (U) , turbulent kinetic energy and streamwise 
turbulent shear stress component (~u’v’) for blowing ratios M = 0.4 
and 0.9. There is a pair of pages for each location. On the first page 
are the run numbers, spanwise positions, and profiles for which the span- 
wise averages were obtained. The dimensional and non-dimensional spanwise- 
averaged profiles are given on the second page. 
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BLOWING REGION AFTER 3 ROWS OF HOLES, X=>148 CM 
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SPANWISE AVERAGED PROFILES M*0.4 START OF RECOVERY AFTER 11 ROWS OF HOLES, X*188 CM 


nOJl-t vO<n 




eg 

O' cc O 3 eg ro 

CsJ 


a 

O O r-~4 eH < . 1 

l/J 

0f'-r-rr)3cog-in3a3'i"0 

\ 


\ 

dCIrgcOgO'-J-g-trcunoO 

> 

OOOOOO 

CNJ 

3 : 

Mcjmn j-g’UMA'tffir.oo 

3 

1 1 1 1 I 1 


-O rgC'-Oog-g-O 
g- nn -om -j-o 

■ — • ■ — l p— I > — i ■ — ( CD CD 


ooooooc 

I t I I I I 


~ OOOOOOOOOOOoO 3: 

> I I I I I I I I I I I I o 

3 


O' 

o 


eg 

tn 3 r-o'f--omgwcoo<foo— i 
N -O <t- — I gj — < O' — • O' rfi -O 3 00 eg 
eg oocaHm^-df'Cioooooo 
y 

— rMcgmro-^irom'^'^Jrg — ioo 

eg 

a 


o 

ii 

eg 

x 

_j 

UJ 

a 


D> 

OC 

UJ 

o 

3 


reimegrorn 

OOOOo 

OQQQfl 
O 3rg enc— 
■O^- OUDgj 

U"\"C| -T 3 — < 


."0 rO rn rO m -H e~ 

oooooooo 

I 

00000000 
.■DrgNg’ag-^HO 
:> -g eg O' eg O' o 
— «cgrg-g~-.— ir^.3 


O' 

o 


II 


ooooooooooooo 


Q 


-g nor — iffiin<n' 30 tD't , go 

_j egoe- O' O' cc l'~-c~r-- coo 
uj OOrg OOOOOOOO—IO 

Q I ■ i ...» t ... I • 

-v. OOOOOOOOOOOOO 


eg 


> 

3 


OcDegegcgcg 

aooooo 

l I I i I I 

aooaoa 

f~-meg.^c- Q\ 
O' CO — . — I — 1 3 


egegegcgrei^ro 

ooooooo 

I I I I I I 

ooooooo 

OOr— cg 301 .nO 
cgrg-H —en —to 


OOOOOOOOOOOOO 

I t I I I I I I I I I I 


l/} 

uj reio^O.^gine-aj-oroaoegcoQ* 

3 LL OOO— I— t— •— I— .OCOgD— IQ 

— ' 3 3—»~. .3—.— . — ■ . — . — «oooo 

u. N 

O O OOOOOOOOOOOOO 

o C 

a. 



o 

< 

ac 

UJ 

:» 

<r 


<t 

z 

o 


UJ 

X 

o 


X 

o 

>- 


00 -O'Oeg p ^>} l CT'C 0 eDao»oegeg 
Oegr^OO'pf t\|-0<ge-0(f|C- X 
O'l'-o^e- egg- coegocj' vt-o o 

............. 

aoo'ooo— ‘egeg j-vnmo -3 -g- 

^ I ^ A ^ r ■ 4 ^ ^ ^ t—A • ^ CD 

rsj 

•r 

O'lnOOOOOOOOOOO II 
— ic-ononooooooo i 

'rmorgne-rnn^noo— < 3 

............. uj 

OO— irgcgeDcnvt’ingD o 


< -J-g-u'iinor-'CT'-OC'— <cgg-o 

z u. -j-a.— .ooconxmcgc-oo 
o 3 iniA-oo-O'Cr-c-coa' jioo 


«/> 3 0000000000033 

z 

UJ 

X 

*- 2 

0 3 a 3 <fTC|egooor-g--j-r->J-e'sj- 

1 uj CT'Ocnccnoocoeg— <ra>oeg 

Z O o— 1 iegeg,- 0 '* - minr~ coa'Hg- 

a s 

Z >- 000000000003 — . 


165 



0.9 BLOWING REGION AFTER 7 ROWS OF HOLES, X=168 CM 
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0.9 START OF RECOVERY AFTER 11 ROrfS OF H01ES» X = 188 CM 
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APPENDIX B 


WIRING DIAGRAMS OF THE 3-D TURBULENT FLOW ANALYZER 


Nomenclatur e 


V 

x 

U 

x 



V 

z 

U 

z 


Linearizer outputs. 

Normalized effective velocities, same as u e ff » u e ff 5 


D , D , D , E , 
x y z x 


X,Y,Z 


X , Y ,Z ,X , 

U U U V 



E E , F , F , F Coefficients in the inverted Jorgensen's 

y z x y z 

matrix M. 


Normalized velocities in wire coordinates. 


Y , Z , X , Y , Z Coefficients of the coordinate transfor- 

v v w w w 

mation matrix N. 


u,v,w 


Normalized velocities in laboratory coordinates. 
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FLOW COMPUTER. MULTIPLIERS AND INVERTERS 
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PLOW COMPUTER SQUARE ROOTS 
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